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INTRODUCTION: 

An  increased  risk  of  pancreatic  cancer  and  reduced  survival  is  seen  among  patients  with  altered 
systemic  metabolism,  including  obesity,  hyperglycemia,  and  hyperinsulinemia.  However,  the 
mechanisms  by  which  altered  metabolism  may  promote  pancreatic  cancer  development  and 
growth  remain  to  be  identified.  The  primary  purpose  of  the  funded  research  is  to  understand  the 
metabolic  alterations  that  promote  development  of  early  pancreatic  cancers  and  impact  their 
proliferation.  To  explore  altered  metabolism  is  patients  with  pancreatic  cancer,  we  are 
conducting  studies  in  participants  from  four  large  prospective  cohorts,  where  individuals  provide 
extensive  data  on  metabolic  phenotypes,  such  as  obesity  and  diabetes,  and  banked  plasma 
samples  for  interrogation.  The  potential  impact  of  understanding  the  mechanisms  underlying 
early  pancreatic  cancer  growth  is  substantial,  with  an  improved  ability  to:  (a)  diagnose 
pancreatic  cancer  at  an  earlier  stage  when  cure  is  possible,  (b)  formulate  preventative 
recommendations  based  on  an  individual’s  metabolic  phenotype,  and  (c)  identify  novel 
treatment  strategies  that  disrupt  pancreatic  tumor  metabolism. 

KEYWORDS: 

Pancreatic  cancer;  Metabolism;  Early  detection 

ACCOMPLISHMENTS: 

What  were  the  major  goals  of  the  project? 

Four  major  tasks  were  described  in  the  approved  Statement  of  Work. 

Major  Task  1 :  For  all  three  Specific  Aims,  local  IRB  and  HRPO  approval  will  be  obtained  to  use 
the  previously  collected  plasma  samples  from  subjects  in  four  large  prospective  cohorts. 

Major  Task  1  included  local  IRB  and  HRPO  approvals  to  be  completed  in  the  first  5  months  of 
the  award  to  allow  evaluation  of  metabolism  in  human  subjects.  This  task  was  completed  in  the 
described  timeframe. 

Major  Task  2:  Identify  appropriate  matched  case-control  set  from  four  cohorts  and  send  plasma 
samples  to  the  Broad  Institute  at  MIT  and  Harvard  University. 

Major  Task  2  included  generation  of  a  human  pancreatic  cancer  case-control  population  and 
aliquotting  and  shipment  of  plasma  samples  for  metabolite  profiling  to  the  Broad  Institute  of  MIT 
and  Harvard  University  in  months  6-9  of  the  award.  This  task  was  completed  in  the  described 
timeframe. 

Major  Task  3:  Generate  and  analyze  metabolomics  dataset. 

Major  Task  3  included  the  analysis  by  LC-MS  of  plasma  samples  for  circulating  metabolites 
(months  10-15),  quality  control  of  the  metabolite  data  (months  16-17),  and  analysis  of  the 
metabolites  related  to  development  of  early  pancreatic  cancer  (months  18-24).  Two  LC-MS 
platforms  were  proposed  for  this  work.  One  method  had  been  completed  in  the  described 
timeframe.  Work  to  profile  metabolites  using  the  second  method  is  ongoing,  as  described  in  the 
previously  submitted  no  cost  extension  request. 
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Major  Task  4:  Generate  and  analyze  metabolites  related  to  branched  chain  amino  acid  (BCAA) 
metabolism. 

Major  Task  4  included  analyses  of  branched  chain  amino  acid  (BCAA)  catabolic  products 
(months  10-15),  quality  control  of  data  (months  16-17),  and  analysis  of  BCAA  catabolic  products 
with  pancreatic  cancer  development  (months  18-23).  Evaluation  of  BCAAs  has  been  completed 
in  the  described  timeframe.  Some  of  the  BCAA  catabolic  products  are  included  on  the  second 
LC-MS  platform  and  their  analysis  is  ongoing. 

What  was  accomplished  under  these  goals? 

1 )  Major  activities:  The  primary  activities  conducted  during  the  second  1 2  months  of  this  award 
were  related  to  analysis  of  circulating  plasma  metabolites  and  their  association  with  pancreatic 
cancer  development  and  progression. 

2)  Specific  objectives:  The  primary  objective  in  the  second  year  of  this  award  was  to  generate 
and  analyze  a  high-quality,  thoroughly-vetted  dataset  of  pancreatic  cancer  cases  and  controls 
with  extensive  clinical  characterization  and  >4,000  measured  plasma  metabolites.  This  dataset 
forms  the  basis  for  defining  metabolic  determinant  of  early  pancreatic  cancer  development  and 
progression. 

3)  Significant  results:  The  important  result  from  the  second  year  of  the  award  is  the  successful 
analysis  of  >2,000  plasma  metabolites  using  our  non-targeted  LC-MS  platforms  in  1,500 
pancreatic  cancer  cases  and  controls,  as  proposed  in  our  award  application. 

Using  LC-MS  Method  1,  we  identified  2,477  circulating  metabolites  that  could  be  reliably  aligned 
across  the  plasma  samples  from  participants  in  four  prospective  cohort  studies  (HPFS,  NHS, 
PHS  and  WHI).  Of  these  metabolites,  1 ,228  were  excluded  because  their  coefficients  of 
variance  were  >  25%.  Of  the  remaining  metabolites,  32  were  excluded  because  data  were 
missing  for  >10%  of  participants.  Finally,  we  assessed  stability  with  processing  delay.  Samples 
from  HPFS  and  NHS  were  shipped  overnight  on  ice.  Thus,  metabolites  not  stable  during  this 
time  period  must  be  excluded  from  analyses.  We  performed  a  pilot  study,  whereby  we  simulated 
these  conditions  and  included  the  samples  in  a  blinded  fashion  within  the  actual  study 
participant  samples.  Of  the  metabolites  that  passed  our  first  two  QC  metrics,  another  1 99 
metabolites  were  not  stable  on  processing  delay  and  were  excluded  from  the  final  data  set. 
Thus,  after  our  3-level  QC  process,  1,018  metabolites  were  of  sufficient  quality  to  be  analyzed  in 
our  nested  case-control  set. 

Using  conditional  logistic  regression,  we  have  identified  20  metabolites  with  FDR  P-value<0.1. 
The  raw  P-values  for  these  metabolites  range  from  3.3  xIO'5  to  1.9  x  10'3.  Of  these  20 
metabolites,  four  are  “knowns”,  including  isoleucine,  leucine,  valine,  and  glycine.  We  identified 
the  three  BCAAs  (isoleucine,  leucine,  and  valine)  in  a  previously  published  study  (Mayers  et  al, 
Nat  Med,  2014).  In  addition  to  these  4  known  metabolites,  we  also  identified  16  “unknowns”, 
which  are  metabolites  of  unknown  composition.  We  are  working  to  determine  the  identity  of 
each  of  these  unknowns.  Interestingly,  seven  of  the  17  non-BCAA  metabolites  have  low 
correlation  coefficients  (rho<0.3)  with  the  BCAAs  or  previously  measured  markers  of  insulin 
resistance,  suggesting  novel  mechanisms  of  association.  As  expected  risk  estimates  for  these 
metabolites  do  not  change  when  considering  known  risk  factors  for  PDAC  in  multivariable- 
adjusted  models,  including  BCAAs,  body-mass  index,  and  circulating  markers  of  insulin 
resistance. 
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In  addition  to  single  metabolite  models,  we  are  currently  building  multi-metabolite  risk  models 
that  simultaneously  incorporate  information  from  correlated  metabolites.  The  goal  of  this  work  is 
to  better  stratify  the  population  by  PDAC  risk,  while  building  a  better  understanding  of  the 
biology  underlying  disease  development. 

4)  Other  achievements:  During  the  second  year  of  this  award,  we  published  four  further 
manuscripts  related  to  the  proposed  work. 

In  the  first  manuscript,  we  examined  the  association  of  plasma  levels  of  25-hydroxyvitamin  D 
(25[OH]D)  with  patient  survival,  demonstrating  that  low  levels  of  25(OH)  are  associated  with 
lower  survival  time.  This  manuscript  was  published  in  Journal  of  Clinical  Oncology. 

Yuan  C,  Qian  ZR,  Babic  A,  Morales-Oyarvide  V,  Rubinson  DA,  Kraft  P,  Ng  K,  Bao  Y, 
Giovannucci  EL,  Ogino  S,  Stampfer  MJ,  Gaziano  JM,  Sesso  HD,  Buring  JE,  Cochrane  BB, 
Chlebowski  RT,  Snetselaar  LG,  Manson  JE,  Fuchs  CS,  Wolpin  BM.  Prediagnostic  Plasma  25- 
Hydroxyvitamin  D  and  Pancreatic  Cancer  Survival.  J  Clin  Oncol.  2016;34(24):2899-905. 

Purpose:  Although  vitamin  D  inhibits  pancreatic  cancer  proliferation  in  laboratory  models,  the 
association  of  plasma  25-hydroxyvitamin  D  (25[OH]D)  with  patient  survival  is  largely 
unexplored.  Patients  and  Methods:  We  analyzed  survival  among  493  patients  from  five 
prospective  U.S.  cohorts  diagnosed  with  pancreatic  cancer  in  1984-2008.  We  estimated  hazard 
ratios  (HRs)  for  death  by  plasma  level  of  25(OH)D  (insufficient,  <20  ng/ml;  relative  insufficiency, 
20-<30  ng/ml;  sufficient  >30  ng/ml)  using  Cox  proportional  hazards  models  adjusted  forage, 
cohort,  race/ethnicity,  smoking,  diagnosis  year,  stage,  and  blood  collection  month.  We  also 
evaluated  33  tagging  single  nucleotide  polymorphisms  (SNPs)  in  the  vitamin  D  receptor  gene 
( VDR ),  requiring  P<0.002  (0.05/33)  for  statistical  significance.  Results:  The  mean  prediagnostic 
plasma  level  of  25(OH)D  was  24.6  ng/ml,  and  165  (33%)  patients  were  vitamin  D  insufficient. 
Compared  to  patients  with  insufficient  levels,  multivariable-adjusted  HRs  for  death  were  0.79 
(95%  Cl,  0.48-1.29)  for  patients  with  relative  insufficiency  and  0.66  (95%  Cl,  0.49-0.90)  for 
patients  with  sufficient  levels  (P-trend=0.01;  see  Table  2  of  attached  manuscript).  These  results 
were  unchanged  after  further  adjustment  for  body-mass  index  and  history  of  diabetes  (P- 
trend=0.02;  see  Table  2  of  attached  manuscript).  The  association  was  strongest  among  patients 
with  blood  collected  within  5  years  of  diagnosis,  with  a  HR  of  0.58  (95%  Cl,  0.35-0.98) 
comparing  patients  with  sufficient  to  patients  with  insufficient  25(OH)D  levels  (see  Table  4  of 
attached  manuscript).  No  SNP  at  VDR  met  our  corrected  significance  threshold  of  P<0.002; 
rs7299460  was  mostly  strongly  associated  with  survival  (HR  per  minor  allele,  0.80;  95%  Cl, 
0.68-0.95;  P=0.01 ;  see  Table  A3  of  attached  manuscript).  Conclusion:  We  observed  longer 
overall  survival  in  patients  with  pancreatic  cancer  who  had  sufficient  prediagnostic  plasma  levels 
of  25(OH)D. 

In  the  second  manuscript,  we  investigated  markers  of  obesity,  leptin  signaling  and  pancreatic 
cancer  risk  in  our  nested  cases  and  controls,  demonstrating  elevated  risk  for  pancreatic  cancer 
among  men  with  high  plasma  leptin  levels.  This  manuscript  was  published  in  Cancer  Research. 

Babic  A,  Bao  Y,  Qian  ZR,  Yuan  C,  Giovannucci  EL,  Aschard  H,  Kraft  P,  Amundadottir  LT, 
Stolzenberg-Solomon  RZ,  Morales-Oyarvide  V,  Ng  K,  Stampfer  MJ,  Ogino  S,  Buring  JE,  Sesso 
HD,  Gaziano  JM,  Rifai  N,  Poliak  MN,  Anderson  ML,  Cochrane  BB,  Luo  J,  Manson  JE,  Fuchs 
CS,  Wolpin  B.  Pancreatic  cancer  risk  associated  with  prediagnostic  plasma  levels  of  leptin  and 
leptin  receptor  genetic  polymorphisms.  Cancer  Res.  2016  Oct  25.  [Epub  ahead  of  print] 
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Leptin  is  an  adipokine  involved  in  regulating  energy  balance  which  has  been  identified  as  a 
potential  biological  link  in  development  of  obesity-associated  cancers  such  as  pancreatic 
cancer.  In  this  prospective,  nested  case-control  study  of  470  cases  and  1094  controls  from  five 
U.S.  cohorts,  we  used  conditional  logistic  regression  to  evaluate  pancreatic  cancer  risk  by 
prediagnostic  plasma  leptin,  adjusting  for  race/ethnicity,  diabetes,  body-mass  index,  physical 
activity,  plasma  C-peptide,  adiponectin,  and  25-hydroxyvitamin  D.  Due  to  known  differences  in 
leptin  levels  by  gender,  analyses  were  conducted  separately  for  men  and  women.  We  also 
evaluated  associations  between  32  tagging  single  nucleotide  polymorphisms  (SNPs)  in  the 
leptin  receptor  (LEPR)  gene  and  pancreatic  cancer  risk.  Leptin  levels  were  higher  in  female 
versus  male  control  participants  (median,  20.8  vs.  6.7pg/mL;  P<0.0001).  Among  men,  plasma 
leptin  was  positively  associated  with  pancreatic  cancer  risk,  and  those  in  the  top  quintile  had  a 
multivariable-adjusted  odds  ratio  (OR)  of  3.02  compared  to  men  in  the  bottom  quintile  (95%  Cl, 
1.27-7.16;  Ptrend=0.02;  see  Table  2  of  attached  manuscript).  Among  women,  circulating  leptin 
was  not  associated  with  pancreatic  cancer  risk  (Ptrend=0.21 ;  see  Table  2  of  attached 
manuscript).  Results  were  similar  across  cohorts  (Pheterogeneity=0.88  for  two  male  cohorts 
and  0.35  for  three  female  cohorts;  see  Figure  2  of  attached  manuscript).  In  genetic  analyses, 
rsl  0493380  in  LEPR  was  associated  with  increased  pancreatic  cancer  risk  among  women,  with 
an  OR  per  minor  allele  of  1 .54  (95%  Cl,  1.1 8-2.02;  multiple  hypothesis-corrected  P=0.03;  see 
Table  3  of  attached  manuscript).  No  SNPs  were  significantly  associated  with  risk  in  men.  In 
conclusion,  higher  prediagnostic  levels  of  plasma  leptin  were  associated  with  an  elevated  risk  of 
pancreatic  cancer  among  men,  but  not  among  women. 

In  the  third  manuscript,  we  worked  with  colleagues  at  Massachusetts  Institute  of  Technology  to 
understand  how  BCAAs  are  used  by  pancreatic  tumors.  This  manuscript  was  published  in 
Science. 

Mayers  JR,  Torrence  ME,  Danai  LV,  Papagiannakopoulos  T,  Davidson  SM,  Bauer  MR,  Lau  AN, 
Ji  BW,  Dixit  PD,  Hosios  AM,  Muir  A,  Chin  CR,  Freinkman  E,  Jacks  T,  Wolpin  BM,  Vitkup  D, 
Vander  Heiden  MG.  Tissue  of  origin  dictates  branched-chain  amino  acid  metabolism  in  mutant 
Kras-driven  cancers.  Science.  2016;353(6304):1 161-5. 

Tumor  genetics  guides  patient  selection  for  many  new  therapies,  and  cell  culture  studies  have 
demonstrated  that  specific  mutations  can  promote  metabolic  phenotypes.  However,  whether 
tissue  context  defines  cancer  dependence  on  specific  metabolic  pathways  is  unknown.  Kras 
activation  and  Trp53  deletion  in  the  pancreas  or  the  lung  result  in  pancreatic  ductal 
adenocarinoma  (PDAC)  or  non-small  cell  lung  carcinoma  (NSCLC),  respectively,  but  despite  the 
same  initiating  events,  these  tumors  use  branched-chain  amino  acids  (BCAAs)  differently  (see 
Figure  1  of  attached  manuscript).  NSCLC  tumors  incorporate  free  BCAAs  into  tissue  protein  and 
use  BCAAs  as  a  nitrogen  source,  whereas  PDAC  tumors  have  decreased  BCAA  uptake.  These 
differences  are  reflected  in  expression  levels  of  BCAA  catabolic  enzymes  in  both  mice  and 
humans  (see  Figure  3  in  the  attached  manuscript).  Loss  of  Bcatl  and  Bcat2,  the  enzymes 
responsible  for  BCAA  use,  impairs  NSCLC  tumor  formation,  but  these  enzymes  are  not  required 
for  PDAC  tumor  formation,  arguing  that  tissue  of  origin  is  an  important  determinant  of  how 
cancers  satisfy  their  metabolic  requirements  (see  Figure  4  of  attached  manuscript). 

In  the  fourth  manuscript,  we  performed  a  study  to  evaluate  factors  at  the  time  of  blood  collection 
that  impact  metabolite  measurements.  These  factors  are  important  to  consider  in  large  human 
studies  that  incorporate  metabolite  profiling.  This  manuscript  was  published  in  Cancer 
Epidemiology  Biomarkers  &  Prevention. 
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Townsend  MK,  Bao  Y,  Poole  EM,  Bertrand  KA,  Kraft  P,  Wolpin  BM,  Clish  CB,  Tworoger  SS. 
Impact  of  Pre-analytic  Blood  Sample  Collection  Factors  on  Metabolomics.  Cancer  Epidemiol 
Biomarkers  Prev.  2016;25(5):823-9. 

BACKGROUND:  Many  epidemiologic  studies  are  using  metabolomics  to  discover  markers  of 
carcinogenesis.  However,  limited  data  are  available  on  the  influence  of  pre-analytic  blood 
collection  factors  on  metabolite  measurement.  METHODS:  We  quantified  166  metabolites  in 
archived  plasma  from  423  Health  Professionals  Follow-up  Study  and  Nurses'  Health  Study 
participants  using  liquid  chromatography-tandem  mass  spectrometry  (LC-MS).  We  compared 
multivariable-adjusted  geometric  mean  metabolite  LC-MS  peak  areas  across  fasting  time, 
season  of  blood  collection,  and  time  of  day  of  blood  collection  categories.  RESULTS:  The 
majority  of  metabolites  (160  of  166  metabolites)  had  geometric  mean  peak  areas  that  were 
within  15%  comparing  samples  donated  after  fasting  9  to  12  versus  >13  hours;  greater 
differences  were  observed  in  samples  donated  after  fasting  <4  hours  (see  Table  3  in  attached 
manuscript).  Metabolite  peak  areas  generally  were  similar  across  season  of  blood  collection, 
although  levels  of  certain  metabolites  (e.g.,  bile  acids  and  purines/pyrimidines)  tended  to  be 
different  in  the  summer  versus  winter  months  (see  Table  4  of  the  attached  manuscript).  After 
adjusting  for  fasting  status,  geometric  mean  peak  areas  for  bile  acids  and  vitamins,  but  not 
other  metabolites,  differed  by  time  of  day  of  blood  collection  (see  Table  5  of  the  attached 
manuscript).  CONCLUSION:  Fasting,  season  of  blood  collection,  and  time  of  day  of  blood 
collection  were  not  important  sources  of  variability  in  measurements  of  most  metabolites  in  our 
study.  However,  considering  blood  collection  variables  in  the  design  or  analysis  of  studies  may 
be  important  for  certain  specific  metabolites,  particularly  bile  acids,  purines/pyrimidines,  and 
vitamins.  IMPACT:  These  results  may  be  useful  for  investigators  formulating  analysis  plans  for 
epidemiologic  metabolomics  studies,  including  determining  which  metabolites  to  a  priori 
exclude  from  analyses. 


What  opportunities  for  training  and  professional  development  has  the  project  provided? 

During  the  second  year  of  this  award,  Dr.  Wolpin  has  met  regularly  with  Dr.  Clary  Clish  at  the 
Broad  Institute  of  MIT  and  Harvard  University  and  Dr.  Peter  Kraft  of  the  Harvard  School  of 
Public  Health.  These  meeting  have  entailed  detailed  review  of  study  design,  analysis  plans,  and 
the  analytic  characteristics  of  the  LC-MS  platforms.  The  close  proximity  of  Dana-Faber  Cancer 
Institute,  Harvard  School  of  Public  Health  and  Broad  Institute  enables  these  regular,  in-person 
meetings.  Dr.  Wolpin  meets  regularly  with  the  Designated  Collaborator  for  the  award,  Dr. 
Charles  Fuchs.  Dr.  Fuchs’  office  is  adjacent  to  Dr.  Wolpin’s  allowing  for  both  regularly 
scheduled  meetings  and  numerous  impromptu  discussions  regarding  design  and  analysis 
considerations.  Furthermore,  Dr.  Wolpin  has  a  graduate  student  and  instructor  in  his  group  that 
accompany  him  to  these  meetings,  allowing  for  further  training  and  professional  development 
opportunities  for  more  junior  investigators.  Dr.  Wolpin  continues  attend  a  number  of  relevant 
conferences  at  Dana-Farber  and  other  nearby  Harvard  institutions.  He  presents  his  work  at  the 
Channing  Laboratory  cohort  study  meetings,  obtaining  real-time  feedback  and  mentoring  from 
senior  faculty.  He  attends  the  DFCI  Seminars  in  Oncology  series  and  had  the  opportunity  to 
present  his  work  at  this  Institute-wide  conference.  He  also  traveled  to  several  conference  and 
delivered  talks  on  his  work,  including  the  American  Association  for  Cancer  Research  (AACR) 
Integrative  Molecular  Epidemiology  Workshop  (Boston,  MA),  Cold  Spring  Harbor  Laboratory 
Workshop  on  Vitamin  D  and  Pancreatic  Cancer  (Cold  Spring  Harbor,  NY),  and  NCI  Cohort 
Consortium  annual  meeting  (Rockville,  MD).  Thus,  Dr.  Wolpin  has  had  ample  opportunities  to 
expand  his  training  and  professional  development  during  the  first  year  of  the  award. 

How  were  the  results  disseminated  to  communities  of  interest? 
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Nothing  to  Report. 

What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

During  the  no-cost  extension  period  of  the  award,  we  will  complete  our  metabolite  analyses, 
including  individual  known  and  unknown  metabolites,  BCAA  catabolic  products,  and  pathway 
analysis  modules  developed  at  the  Broad  Institute. 

IMPACT: 

What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

The  initiated  work  will  advance  our  understanding  of  how  pancreatic  cancers  use  nutrients  to 
promote  their  growth.  This  has  important  applications  to  early  detection  of  pancreatic  cancer 
and  identifying  new  treatments  that  cut  off  the  nutrient  supply  to  tumors. 

What  was  the  impact  on  other  disciplines? 

The  work  performed  within  this  award  will  advance  the  field  technically  by  demonstrating  how  to 
pursue  large-scale  metabolism  studies  using  banked  blood  samples  from  thousands  of 
individuals.  This  has  relevance  to  many  disease  states,  beyond  cancer. 

What  was  the  impact  on  technology  transfer? 

Nothing  to  Report. 

What  was  the  impact  on  society  beyond  science  and  technology? 

This  work  has  the  potential  to  increase  awareness  around  pancreatic  cancer,  which  is  now  the 
third  leading  cause  of  cancer-related  death  in  the  United  States.  Furthermore,  this  work  has  the 
potential  to  highlight  the  link  between  obesity,  diabetes,  and  cancer,  which  has  important 
implications  for  lifestyle  recommendations  to  the  general  public. 

CHANGES/PROBLEMS: 

Changes  in  approach  and  reasons  for  change 

Nothing  to  Report 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Nothing  to  Report 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  Report 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards, 
and/or  select  agents 

Nothing  to  Report 

Significant  changes  in  use  or  care  of  human  subjects 

Nothing  to  Report 

Significant  changes  in  use  or  care  of  vertebrate  animals. 

Nothing  to  Report 

Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  Report 
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PRODUCTS: 


Publications,  conference  papers,  and  presentations 

Journal  publications. 

1 .  Yuan  C,  Qian  ZR,  Babic  A,  Morales-Oyarvide  V,  Rubinson  DA,  Kraft  P,  Ng  K,  Bao  Y, 
Giovannucci  EL,  Ogino  S,  Stampfer  MJ,  Gaziano  JM,  Sesso  HD,  Buring  JE,  Cochrane  BB, 
Chlebowski  RT,  Snetselaar  LG,  Manson  JE,  Fuchs  CS,  Wolpin  BM.  Prediagnostic  Plasma 
25-Hydroxyvitamin  D  and  Pancreatic  Cancer  Survival.  J  Clin  Oncol.  2016;34(24):2899-905. 

2.  Babic  A,  Bao  Y,  Qian  ZR,  Yuan  C,  Giovannucci  EL,  Aschard  H,  Kraft  P,  Amundadottir  LT, 
Stolzenberg-Solomon  RZ,  Morales-Oyarvide  V,  Ng  K,  Stampfer  MJ,  Ogino  S,  Buring  JE, 
Sesso  HD,  Gaziano  JM,  Rifai  N,  Poliak  MN,  Anderson  ML,  Cochrane  BB,  Luo  J,  Manson 
JE,  Fuchs  CS,  Wolpin  BM.  Pancreatic  cancer  risk  associated  with  prediagnostic  plasma 
levels  of  leptin  and  leptin  receptor  genetic  polymorphisms.  Cancer  Res.  2016  Oct  25.  [Epub 
ahead  of  print] 

3.  Mayers  JR,  Torrence  ME,  Danai  LV,  Papagiannakopoulos  T,  Davidson  SM,  Bauer  MR,  Lau 
AN,  Ji  BW,  Dixit  PD,  Hosios  AM,  Muir  A,  Chin  CR,  Freinkman  E,  Jacks  T,  Wolpin  BM, 
Vitkup  D,  Vander  Heiden  MG.  Tissue  of  origin  dictates  branched-chain  amino  acid 
metabolism  in  mutant  Kras-driven  cancers.  Science.  2016;353(6304):1 161-5. 

4.  Townsend  MK,  Bao  Y,  Poole  EM,  Bertrand  KA,  Kraft  P,  Wolpin  BM.  Clish  CB,  Tworoger 
SS.  Impact  of  Pre-analytic  Blood  Sample  Collection  Factors  on  Metabolomics.  Cancer 
Epidemiol  Biomarkers  Prev.  2016;25(5):823-9 

Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  Report. 

Other  publications,  conference  papers,  and  presentations. 

Local  presentations: 


2016  Precision  Medicine  for  Pancreatic 

Adenocarcinoma 

Dana-Farber  Cancer  Institute 


Center  for  Cancer 
Precision  Medicine  Tumor 
Board  Conference 


2016  Precision  Medicine  for  Pancreatic 

Adenocarcinoma 

Dana-Farber  Cancer  Institute 


Conference  of  the  Hale 
Center  for  Pancreatic 
Cancer  Research 


National  presentations: 

2015  Epistatic  Targeting  of  Primary  and  Metastatic  Invited  Speaker 

Pancreatic  Cancer 

Stand  Up  to  Cancer  -  Cancer  Research  UK  -  Lustgarten  Philadelphia,  PA 
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Foundation  Dream  Team  Finalist  Presentation 


2015  Research  Biopsy  Program  in  Patients  with 

Metastatic  Pancreatic  Cancer 

Lustgarten  Foundation  Scientific  Meeting 


2015  The  Need  for  Both  Edges  of  the  Sword:  Early 
Detection  and  New  Treatment  Approaches  for 
Pancreatic  Cancer  Mortality  Reduction 

MD  Anderson  Cancer  Center  Ahmad  Center  Seminar  Series 

2016  The  Double-edged  Sword  for  Pancreatic  Cancer 
Mortality  Reduction:  Identifying  New  Approaches 
for  Early  Detection  and  Treatment 

Oregon  Health  &  Science  University  Seminar  Series,  Brenden- 

Colson  Center  for  Pancreatic  Care 

2016  Altered  Systemic  Metabolism  and  Pancreatic 

Adenocarcinoma 

University  of  California  at  San  Francisco,  Pancreas  Center  and 

Diabetes  &  Obesity  Research  Seminar  Series 

2016  Pancreatic  Adenocarcinoma  and  Altered  Host 

Metabolism 

American  Society  of  Preventative  Oncology  40th  Annual 

Meeting 

2016  Dana-Farber  Pancreatic  Cancer  Research  Tissue 

Platform 

Pancreatic  Cancer  Action  Network  RAN-2  Grant  Meeting 

2016  Defining  Novel  Approaches  to  Early  Detection 

and  Treatment  of  Pancreatic  Adenocarcinoma 

The  Ohio  State  University  Comprehensive  Cancer  Center 

Seminar  Series 

2016  Defining  Novel  Approaches  to  Early  Detection 

and  Treatment  of  Pancreatic  Adenocarcinoma 

University  of  Michigan  Comprehensive  Cancer  Center  Seminar 

Series 

2016  Metabolomics  and  Cancer 


Invited  Speaker 

Cold  Spring  Harbor, 
NY 

Invited  Speaker 

Houston,  TX 
Invited  Speaker 

Portland,  OR 

Invited  Speaker 
San  Francisco,  CA 

Invited  Speaker 
Columbus,  OH 

Invited  Speaker 
Rochester,  NY 

Invited  Speaker 
Columbus,  OH 

Invited  Speaker 
Ann  Arbor,  Ml 

Invited  Speaker 
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American  Association  for  Cancer  Research  (AACR)  Integrative  Boston,  MA 
Molecular  Epidemiology  Workshop 

International  Invited  Presentations 

2016  Altered  Systemic  Metabolism  and  Early  Invited  Speaker 

Pancreatic  Adenocarcinoma 

Pezcoller  Symposium  Trento,  Italy 

In  the  above  presentations,  metabolism  changes  due  to  pancreatic  cancer  were  discussed, 
including  how  these  changes  can  inform  early  detection  and  prognosis  of  pancreatic  cancer. 

Website(s)  or  other  Internet  site(s) 

Nothing  to  Report. 

Technologies  or  techniques 

Nothing  to  Report. 

inventions,  patent  applications,  and/or  licenses 

Nothing  to  Report. 

Other  Products 

Nothing  to  Report. 

PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 

What  individuals  have  worked  on  the  project? 

Brian  Wolpin.  No  change. 

Chen  Yuan.  No  change. 

Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key 
personnel  since  the  last  reporting  period? 

See  attached  Other  Support  page. 

What  other  organizations  were  involved  as  partners? 

Nothing  to  Report. 

SPECIAL  REPORTING  REQUIREMENTS 

COLLABORATIVE  AWARDS: 

Not  applicable. 

QUAD  CHARTS: 

Not  applicable. 

APPENDICES: 

Four  manuscripts: 
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1 .  Yuan  C,  Qian  ZR,  Babic  A,  Morales-Oyarvide  V,  Rubinson  DA,  Kraft  P,  Ng  K,  Bao  Y, 
Giovannucci  EL,  Ogino  S,  Stampfer  MJ,  Gaziano  JM,  Sesso  HD,  Buring  JE,  Cochrane  BB, 
Chlebowski  RT,  Snetselaar  LG,  Manson  JE,  Fuchs  CS,  Wolpin  BM.  Prediagnostic  Plasma 
25-Hydroxyvitamin  D  and  Pancreatic  Cancer  Survival.  J  Clin  Oncol.  2016;34(24):2899-905. 

2.  Babic  A,  Bao  Y,  Qian  ZR,  Yuan  C,  Giovannucci  EL,  Aschard  H,  Kraft  P,  Amundadottir  LT, 
Stolzenberg-Solomon  RZ,  Morales-Oyarvide  V,  Ng  K,  Stampfer  MJ,  Ogino  S,  Buring  JE, 
Sesso  HD,  Gaziano  JM,  Rifai  N,  Poliak  MN,  Anderson  ML,  Cochrane  BB,  Luo  J,  Manson 
JE,  Fuchs  CS,  Wolpin  BM.  Pancreatic  cancer  risk  associated  with  prediagnostic  plasma 
levels  of  leptin  and  leptin  receptor  genetic  polymorphisms.  Cancer  Res.  2016  Oct  25.  [Epub 
ahead  of  print] 

3.  Mayers  JR,  Torrence  ME,  Danai  LV,  Papagiannakopoulos  T,  Davidson  SM,  Bauer  MR,  Lau 
AN,  Ji  BW,  Dixit  PD,  Hosios  AM,  Muir  A,  Chin  CR,  Freinkman  E,  Jacks  T,  Wolpin  BM, 
Vitkup  D,  Vander  Heiden  MG.  Tissue  of  origin  dictates  branched-chain  amino  acid 
metabolism  in  mutant  Kras-driven  cancers.  Science.  2016;353(6304):1 161-5. 

4.  Townsend  MK,  Bao  Y,  Poole  EM,  Bertrand  KA,  Kraft  P,  Wolpin  BM.  Clish  CB,  Tworoger 
SS.  Impact  of  Pre-analytic  Blood  Sample  Collection  Factors  on  Metabolomics.  Cancer 
Epidemiol  Biomarkers  Prev.  2016;25(5):823-9 
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Purpose 

Although  vitamin  D  inhibits  pancreatic  cancer  proliferation  in  laboratory  models,  the  association  of 
plasma  25-hydroxyvitamin  D  [25(OH)D]  with  patient  survival  is  largely  unexplored. 

Patients  and  Methods 

We  analyzed  survival  among  493  patients  from  five  prospective  US  cohorts  who  were  diagnosed 
with  pancreatic  cancer  from  1984  to  2008.  We  estimated  hazard  ratios  (HRs)  for  death  by  plasma 
level  of  25(OH)D  (insufficient,  <  20  ng/mL;  relative  insufficiency,  20  to  <  30  ng/mL;  sufficient  a  30 
ng/mL)  by  using  Cox  proportional  hazards  regression  models  adjusted  for  age,  cohort,  race  and 
ethnicity,  smoking,  diagnosis  year,  stage,  and  blood  collection  month.  We  also  evaluated  30  tagging 
single-nucleotide  polymorphisms  in  the  vitamin  D  receptor  gene,  requiring  P<  .002  (0.05  divided  by 
30  genotyped  variants)  for  statistical  significance. 

Results 

Mean  prediagnostic  plasma  level  of  25(OH)D  was  24.6  ng/mL,  and  1 65  patients  (33%)  were  vitamin 
D  insufficient.  Compared  with  patients  with  insufficient  levels,  multivariable-adjusted  HRs  for  death 
were  0.79  (95%  Cl,  0.48  to  1.29)  for  patients  with  relative  insufficiency  and  0.66  (95%  Cl,  0.49  to 
0.90)  for  patients  with  sufficient  levels  (Ptrend  =  .01).  These  results  were  unchanged  after  further 
adjustment  for  body  mass  index  and  history  of  diabetes  ( P  trend  =  .02).  The  association  was 
strongest  among  patients  with  blood  collected  within  5  years  of  diagnosis,  with  an  HR  of  0.58  (95% 
Cl,  0.35  to  0.98)  comparing  patients  with  sufficient  to  patients  with  insufficient  25(OH)D  levels.  No 
single-nucleotide  polymorphism  at  the  vitamin  D  receptor  gene  met  our  corrected  significance 
threshold  of  P <  .002;  rs7299460  was  most  strongly  associated  with  survival  (HR  per  minor  allele, 
0.80;  95%  Cl,  0.68  to  0.95;  P=  .01). 

Conclusion 

We  observed  longer  overall  survival  in  patients  with  pancreatic  cancer  who  had  sufficient  pre¬ 
diagnostic  plasma  levels  of  25(OH)D. 
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INTRODUCTION 


Pancreatic  cancer  is  the  third-leading  cause  of 
cancer-related  deaths  in  the  United  States,  and 
most  patients  die  within  12  months  of  diagnosis.1 
Other  than  disease  stage  at  presentation,  few 
patient  characteristics  or  circulating  markers  have 
been  identified  that  predict  survival  time  in  pa¬ 
tients  with  pancreatic  ductal  adenocarcinoma.2 

Vitamin  D  is  a  fat-soluble  vitamin  that  can  be 
ingested  in  foods  and  is  produced  endogenously 
when  sunlight  triggers  synthesis  in  the  skin. 
Metabolites  of  vitamin  D  bind  to  the  vitamin  D 


receptor  (VDR)  within  cells,  which  regulate  tran¬ 
scription  of  target  genes  by  interacting  with  vi¬ 
tamin  D  response  elements.’  In  addition  to 
the  essential  roles  of  vitamin  D  in  calcium  ho¬ 
meostasis  and  bone  health,  emerging  evidence 
indicates  that  vitamin  D  and  its  analogs  may 
inhibit  cancer  development  and  growth.3"5  Recent 
studies  in  laboratory  models  of  pancreatic  cancer 
have  demonstrated  the  therapeutic  potential  of 
vitamin  D  analogs  by  engagement  of  VDR  in  tu¬ 
mor  cells  and  supportive  cells  within  the  tumor 
stroma.6"8  Although  studies  that  have  evaluated 
circulating  vitamin  D  and  pancreatic  cancer  in¬ 
cidence  have  been  conflicting,9’10  few  studies  have 
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been  conducted  to  assess  the  association  of  vitamin  D  levels  with 
survival  in  patients  with  pancreatic  cancer.11,12 

To  assess  the  potential  prognostic  effect  of  vitamin  D  status, 
we  evaluated  the  association  between  prediagnostic  plasma  25- 
hydroxyvitamin  D  [25(OH)D]  and  overall  survival  among  patients 
with  pancreatic  cancer  drawn  from  five  large  US  prospective  co¬ 
horts.  In  addition,  we  examined  survival  in  relation  to  common 
germline  genetic  variants  in  VDR. 


PATIENTS  AND  METHODS 


Study  Population 

Our  study  population  included  patients  with  pancreatic  cancer  from 
five  US  prospective  cohort  studies.  The  Health  Professionals  Follow-up 
Study  (HPFS)  was  established  in  1986  when  51,529  male  health  pro¬ 
fessionals,  age  40  to  75  years,  responded  to  a  mailed  questionnaire.  The 
Nurses’  Health  Study  (NHS)  was  initiated  in  1976  when  121,700  female 
registered  nurses,  age  30  to  55  years,  responded  to  a  mailed  questionnaire. 
The  Physicians’  Health  Study  I  (PHS  I)  is  a  randomized  clinical  trial  of 
aspirin  and  (3-carotene  that  enrolled  22,071  male  physicians,  age  40  to 
84  years,  in  1982.  After  completion  of  the  randomly  assigned  components, 
study  participants  were  observed.  The  Women’s  Health  Initiative  (WHI) 
Observational  Study  enrolled  93,676  postmenopausal  women,  age  50  to 
79  years,  between  1994  and  1998.  The  Women’s  Health  Study  (WHS)  is 
a  randomized,  clinical  trial  of  low-dose  aspirin  and  vitamin  E  that  enrolled 
39,876  female  professionals,  age  a  45  years,  between  1992  and  1995.  The 
trial  was  completed  in  2004,  and  participants  were  observed.  Details  of 
these  cohorts  have  been  described  previously  (HPFS,13  NHS,14  PHS,15 
WHI,16  and  WHS17).  The  study  was  approved  by  the  Human  Research 
Committee  at  the  Brigham  and  Women’s  Hospital,  and  participants 
provided  informed  consent. 

We  included  493  patients  with  pancreatic  adenocarcinoma  who  were 
diagnosed  through  2008  and  who  had  available  plasma  and  no  prior 
history  of  cancer,  except  nonmelanoma  skin  cancer.  Incident  cases  of 
pancreatic  cancer  were  identified  by  self-report  or  during  follow-up  of 
participant  deaths.  Deaths  were  ascertained  from  next-of-kin  or  the  US 
Postal  Service  and  by  searching  the  National  Death  Index;  this  method  has 
been  shown  to  capture  >  98%  of  deaths.18  Diagnoses  were  confirmed  by 
review  of  medical  records,  death  certificates,  and/or  tumor  registry  data  by 
study  physicians  who  were  blinded  to  exposure  data.  Patients  with  non¬ 
adenocarcinoma  histology  or  unclear  survival  time  were  excluded. 

Individual  characteristics  and  habits  were  obtained  from  baseline 
questionnaires  at  study  enrollment  in  PHS,  WHI,  and  WHS,  and  from 
questionnaires  that  preceded  the  date  of  blood  draw  in  HPFS  and  NHS.  In 
all  cohorts,  data  were  available  for  age  at  blood  draw,  sex,  race/ethnicity, 
weight,  height,  smoking  status,  and  history  of  diabetes.  Date  of  diagnosis 
and  pancreatic  cancer  stage  at  diagnosis  were  obtained  from  medical  record 
review.  Cancer  stage  was  classified  as  local  disease  amenable  to  surgical 
resection;  locally  advanced  disease  with  extrapancreatic  extension  ren¬ 
dering  it  unresectable,  but  without  distant  metastases;  distant  metastatic 
disease;  or  unknown. 

Plasma  Samples 

Blood  samples  were  collected  from  18,225  male  patients  in  HPFS 
from  1993  to  1995;  32,826  female  patients  in  NHS  from  1989  to  1990; 
14,916  male  patients  in  PHS  from  1982  to  1984;  93,676  female  patients  in 
WHI  from  1994  to  1998;  and  28,345  female  patients  in  WHS  from  1992  to 
1995.  All  blood  samples  were  continuously  stored  in  well-monitored 
freezers.  Details  on  blood  draw  procedures,  transportation,  and  storage 
of  plasma  samples  in  these  cohorts  have  been  described  elsewhere.19  As 
described  previously,  plasma  levels  of  25(OH)D  were  assayed  in  the 
laboratory  of  Nader  Rifai,  MD  (Children’s  Hospital,  Boston,  MA),  by 
using  the  25-Hydroxyvitamin  D  Enzyme  Immunoassay  Kit  from 


Immunodiagnostic  Systems  (Tyne  8c  Wear,  United  Kingdom)  per  man¬ 
ufacturer  instructions.9  Mean  intra-assay  coefficients  of  variance  were  < 
9%  for  blinded,  replicate,  quality  control  samples.9 

Single  Nucleotide  Polymorphism  Selection  and  Genotyping 

We  selected  36  tagging  single  nucleotide  polymorphisms  (SNPs)  in 
the  VDR  gene  ±  20  kb  by  using  the  tagger  algorithm  implemented  in 
Haploview,  with  a  cutoff  of  r2  =  0.8  and  a  minor  allele  frequency  of  a  5% 
in  Whites  from  the  HapMap  Project  database,  and  we  forced  in  one 
previously  reported  SNP  related  to  prostate  cancer  risk  (rsl  15  74  1  43).20  Of 
397  patients  with  pancreatic  cancer,  DNA  was  extracted  from  archived 
huffy  coat  samples  with  Qiagen  QIAmp  (Valencia,  CA)  and  whole-genome 
amplified  with  GE  Healthcare  Genomiphi  (Pittsburgh,  PA).  All  genotyping 
was  carried  out  at  the  Partners  Healthcare  Center  for  Personalized  Genetic 
Medicine  by  using  a  custom-designed  Illumina  Golden  Gate  genotyping 
assay  (San  Diego,  CA).  Three  tagging  SNPs  were  not  supported  by  the 
Golden  Gate  platform  and  so  could  not  be  genotyped.  Three  SNPs  deviated 
from  Hardy- Weinberg  Equilibrium  at  P  <  .01  and  were  excluded.  Replicate 
samples  included  for  quality  control  (n  =  44  sample  groups)  had  mean 
genotype  concordance  of  97.8%  across  the  30  SNPs. 

Statistical  Analysis 

Association  of  plasma  25(OH)D  with  overall  survival  was  examined 
by  using  Cox  proportional  hazards  regression  models  to  calculate  hazard 
ratios  (HRs)  and  95%  CL  Overall  survival  time  was  calculated  from  the 
date  of  cancer  diagnosis  until  the  date  of  death  or  last  follow-up  if 
a  participant  was  still  alive.  Although  optimal  levels  of  25(OH)D  have  not 
been  definitively  determined,  plasma  25(OH)D  <  20  ng/mL  has  been 
defined  as  insufficiency,  20  to  <  30  ng/mL  as  relative  insufficiency,  and 
&  30  ng/mL  as  sufficient21,22;  therefore,  we  investigated  patient  survival  by 
these  categories  of  plasma  25(OH)D.  To  compare  more  extreme  levels  of 
25(OH)D  in  secondary  analyses,  we  assessed  the  association  between 
quintiles  of  25(OH)D  and  survival.  We  examined  HRs  for  each  gender- 
specific  cohort  individually  and  computed  a  summary  HR  by  using  the 
DerSimonian  and  Laird  random  effects  model.23  Heterogeneity  was  tested 
by  using  Cochran’s  Q  statistic.24  Two-sided  tests  for  trend  were  calculated 
by  entering  log-transformed  25(OH)D  as  a  continuous  variable  into  Cox 
proportional  hazards  regression  models,  given  that  25(OH)D  was  not 
normally  distributed.  In  multivariable  models,  we  adjusted  for  potential 
confounders,  including  age  at  diagnosis,  race  and  ethnicity,  smoking  status, 
year  of  diagnosis,  cancer  stage,  month  of  blood  collection,  and  time  be¬ 
tween  blood  collection  and  cancer  diagnosis.  Subsequently,  we  adjusted 
for  body  mass  index  (BMI)  and  history  of  diabetes,  which  were  prog¬ 
nostic  factors  identified  in  prior  studies  of  these  populations.25,26  The 
proportionality  of  hazards  assumption  was  satisfied  by  evaluating  a  time- 
dependent  variable,  which  was  the  cross-product  of  log-transformed 
25(OH)D  and  time  ( P  =  .58). 

We  estimated  median  survival  time  and  survival  curves  for  partic¬ 
ipants  in  each  category  adjusted  for  covariates  by  using  direct  adjusted 
survival  estimation.27,28  This  method  uses  Cox  proportional  hazards  re¬ 
gression  to  estimate  probabilities  of  survival  at  each  follow-up  time  point 
for  each  individual  and  averages  them  to  obtain  an  overall  survival  esti¬ 
mate.  We  assessed  statistical  interaction  by  BMI,  tobacco  use,  cancer  stage, 
season  of  blood  collection,  and  time  between  blood  collection  and  cancer 
diagnosis  by  entering  main  effect  terms  and  a  cross-product  term  of  log- 
transformed  25(OH)D  and  stratification  variable  into  the  model,  and 
evaluated  likelihood  ratio  tests. 

Pancreatic  cancer  and  stromal  cells  express  VDR,  which  translocates 
to  the  nucleus  and  binds  vitamin  D  response  elements  to  regulate  gene 
expression  and  mediate  the  activity  of  vitamin  D.  Thus,  we  examined  the 
association  of  tagging  SNPs  in  VDR  with  survival  by  including  each  three- 
level  genotype  as  a  continuous  variable  (additive  model)  in  multivariable- 
adjusted  Cox  proportional  hazards  regression  models.  SNPs  were 
considered  statistically  significant  if  P  <  .002  (0.05  divided  by  30  genotyped 
variants).  We  subsequently  examined  the  joint  association  of  25(OH)D 
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Vitamin  D  and  Pancreatic  Cancer  Survival 


and  VDR  genotype  in  Cox  regression  models.  All  analyses  were  performed 
with  SAS  9.3  statistical  package.  All  P  values  were  two  sided. 


RESULTS 


Baseline  characteristics  of  patients  with  pancreatic  cancer  by 
category  of  prediagnostic  25(OH)D  and  by  cohort  are  listed  in 
Table  1  and  Appendix  Table  A1  (online  only),  respectively.  Median 
time  between  blood  draw  and  pancreatic  cancer  diagnosis  was 
6.7  years.  Mean  prediagnostic  25(OH)D  was  24.6  ng/mL,  with  33% 
of  patients  classified  as  vitamin  D  insufficient.  Among  those  with 
known  disease  stage,  16.0%  had  localized  disease,  29.6%  had  lo¬ 
cally  advanced  disease,  and  54.3%  had  metastatic  disease.  Median 
survival  by  cancer  stage  was  17  months  for  those  with  localized 
disease,  10  months  for  those  with  locally  advanced  disease,  and 
3  months  for  those  with  metastatic  disease.  At  the  end  of  follow-up, 
464  patients  (94%)  were  deceased. 

Higher  plasma  25(OH)D  levels  were  associated  with  greater 
survival  (P  trend  =  .01;  Table  2  and  Fig  1).  Compared  with  patients 
with  insufficient  levels  of  vitamin  D,  the  multivariable-adjusted 
HRs  for  death  were  0.79  (95%  Cl,  0.48  to  1.29)  for  patients  with 
relative  insufficiency  and  0.66  (95%  Cl,  0.49  to  0.90)  for  patients 
with  sufficient  25(OH)D  levels.  After  adjustment  for  time  between 


Table  1.  Baseline  Characteristics  of  Patients  With  Pancreatic  Cancer  by  Pre¬ 
diagnostic  Plasma  25(OH)D  Level 


Plasma  25(OH)D 

Characteristic 

Insufficient 

Relative  Insufficiency 

Sufficient 

25(OH)D  level 

ng/mL 

<  20 

20  to  <  30 

IV 

CO 

o 

nM 

<  50 

50  to  <  75 

a  75 

No.  of  patients 

165 

212 

116 

Age  at  blood  draw,  years 

63.9  (7.7) 

64.0  (8.8) 

61.6  (8.8) 

Age  at  diagnosis,  years 

71.1  (8.1) 

71.7  (8.1) 

70.2  (8.7) 

Female  sex 

126  (76.4) 

152  (71.7) 

67  (57.8) 

Race/ethnicity 

White 

142  (86.1) 

188  (88.7) 

106  (91.4) 

Black 

14  (8.5) 

4  (1 .9) 

0  (0) 

Other 

5  (3.0) 

8  (3.8) 

2  (1.7) 

Missing 

4  (2.4) 

12  (5.7) 

8  (6.9) 

Body  mass  index,  kg/m2 

27.2  (5.0) 

26.7  (5.5) 

24.9  (3.5) 

Physical  activity,  MET  h/wk 

11.5  (14.7) 

16.8  (20.6) 

26.9  (29.4) 

Diabetes  at  blood  draw 

10  (6.1) 

12  (5.7) 

7  (6.0) 

Tobacco  use 

Never 

65  (39.4) 

96  (45.3) 

43  (37.1) 

Past 

69  (41 .8) 

90  (42.5) 

62  (53.4) 

Current 

29  (17.6) 

25  (1 1 .8) 

11  (9.5) 

Missing 

2  (1 .2) 

1  (0.5) 

0  (0) 

Multivitamin  use 

60  (36.4) 

104  (49.1) 

53  (45.7) 

Alcohol,  >  1  drink/d 

46  (27.9) 

45  (21 .2) 

32  (27.6) 

Diagnosis  period 

1984-2000 

94  (57.0) 

111  (52.4) 

65  (56.0) 

2001-2008 

71  (43.0) 

101  (47.6) 

51  (44.0) 

Cancer  stage 

Localized 

21  (12.7) 

25  (1 1 .8) 

19  (16.4) 

Locally  advanced 

48  (29.1) 

49  (23.1) 

23  (19.8) 

Metastatic 

68  (41 .2) 

97  (45.8) 

55  (47.4) 

Unknown 

28  (17.0) 

41  (19.3) 

19  (16.4) 

NOTE.  Continuous  variables  reported  as  mean  (standard  deviation)  and  cate¬ 
gorical  variables  reported  as  No.  (%),  unless  otherwise  noted. 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  MET,  metabolic  equivalent. 


blood  collection  and  diagnosis,  BMI,  and  history  of  diabetes,  our 
results  remained  largely  unchanged  (Table  2).  In  analyses  to 
compare  more  extreme  values  of  25(OH)D,  a  similar  inverse  re¬ 
lationship  was  noted  between  plasma  25(OH)D  and  patient  sur¬ 
vival  (Table  3).  Compared  with  those  in  the  bottom  quintile  of 
25(OH)D,  patients  in  the  top  quintile  had  a  multivariable-adjusted 
HR  for  death  of  0.68  (95%  Cl,  0.48  to  0.95).  We  considered  the 
possible  influence  of  subclinical  malignancy  on  plasma  25(OH)D. 
After  excluding  19  patients  who  were  diagnosed  with  pancreatic 
cancer  within  one  year  of  blood  collection,  our  results  were  not 
materially  altered  (Appendix  Table  A2,  online  only). 

No  statistically  significant  interactions  were  observed  by  BMI, 
tobacco  use,  cancer  stage,  season  of  blood  collection,  or  time 
between  blood  collection  and  diagnosis  (Table  4).  A  stronger 
association  of  25(OH)D  and  patient  survival  was  identified  in 
patients  with  blood  collected  within  5  years  of  diagnosis.  Among 
patients  with  blood  collected  <  5  years  before  cancer  diagnosis,  HR 
for  death  was  0.58  (95%  Cl,  0.35  to  0.98),  a  comparison  of  those 
with  sufficient  levels  and  those  with  insufficient  levels  of  25(OH)D. 
We  evaluated  the  association  of  25(OH)D  and  patient  survival 
across  the  cohort  study  populations  (Fig  2),  and  noted  no  sta¬ 
tistically  significant  heterogeneity  ( P  heterogeneity  =  .78).  In 
a  comparison  of  patients  with  sufficient  25(OH)D  levels  and 
patients  with  insufficient  levels  of  plasma  vitamin  D,  multivariable- 
adjusted  HRs  for  death  in  sex-specific  meta-analyses  were  0.61 
(95%  Cl,  0.42  to  0.89)  for  females  (NHS,  WHI,  and  WHS)  and 
0.79  (95%  Cl,  0.45  to  1.39)  for  males  (HPFS  and  PHS). 

We  next  examined  genotypic  variation  at  VDR  in  relation  to 
patient  survival.  Three  tagging  SNPs  in  VDR  were  associated  with 
survival  to  P  <  .05  in  an  additive  model  of  inheritance,  but  no  SNP 
met  the  predefined  significance  threshold  after  multiple  hypothe¬ 
sis  testing  correction  (Appendix  Table  A3,  online  only).  The  most 
significant  SNP  was  rs7299460,  which  had  a  per-allele  HR  for  death  of 
0.80  (95%  Cl,  0.68  to  0.95;  P  =  .01).  In  an  exploratory  analysis, 
patients  who  were  homozygous  recessive  at  rs7299460  and  with 
sufficient  25(OH)D  had  the  best  overall  survival,  with  an  HR  for  death 
of  0.45  (95%  Cl,  0.22  to  0.93;  Appendix  Table  A4,  online  only). 


DISCUSSION 


In  this  prospective  study  of  patients  with  pancreatic  cancer  from 
five  large  US  cohorts,  patients  with  sufficient  levels  of  prediagnostic 
25(OH)D  had  a  35%  lower  hazard  for  death  compared  with  those 
who  were  vitamin  D  deficient.  Association  of  prediagnostic 
25(OH)D  with  survival  remained  unchanged  after  consideration  of 
multiple  potential  confounding  factors  and  seemed  to  be  strongest 
when  blood  was  collected  within  five  years  before  pancreatic  cancer 
diagnosis.  Although  several  polymorphisms  at  the  VDR  locus 
were  nominally  associated  with  survival,  no  SNP  met  the  threshold 
for  statistical  significance  after  multiple  hypothesis  testing  correc¬ 
tion.  In  aggregate,  these  results  indicate  that  prediagnostic  plasma 
25(OH)D  is  a  prognostic  factor  in  patients  with  pancreatic  cancer. 

In  laboratory  models,  vitamin  D  and  its  analogs  have  been 
shown  to  induce  differentiation,  promote  apoptosis,  and  in¬ 
hibit  proliferation  of  pancreatic  cancers.3,5,29"31  Furthermore, 
recent  studies  have  indicated  direct  effects  of  ligand  binding  to 
VDR  in  pancreatic  cancer  cells,  but  also  indirect  effects  by  VDR 
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Table  2.  HRs  for  Death  Among  Patients  With  Pancreatic  Cancer  by  Plasma  25(OH)D  Level 

Model 

Plasma  25(OH)D 

/"trend 

Insufficient 

Relative  Insufficiency 

Sufficient 

25(OH)D  level 

ng/mL 

<  20 

20  to  <  30 

o 

CO 

Al 

nM 

<  50 

50  to  <  75 

-  75 

Median  survival,  months 

5 

7 

8 

Person-months 

1,712 

2,843 

1,348 

Patients/deaths 

165/154 

212/1 99 

116/111 

HR  (95%  Cl)t 

1.0 

0.79  (0.48  to  1.29) 

0.66  (0.49  to  0.90) 

.01 

HR  (95%  Cl)t 

1.0 

0.83  (0.52  to  1.32) 

0.67  (0.49  to  0.92) 

.02 

HR  (95%  CDS 

1.0 

0.77  (0.46  to  1.29) 

0.62  (0.44  to  0.86) 

.02 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  HR,  hazard  ratio. 

‘Test  for  trend  calculated  by  entering  log-transformed  25(OH)D  as  a  continuous  variable  in  Cox  proportional  hazards  regression  models. 

tMeta-analysis  of  cohort-specific  HRs  (95%  Cl)  from  Cox  proportional  hazards  regression  models  adjusted  for  age  at  diagnosis,  race  and  ethnicity  (white,  black,  other, 
unknown),  smoking  status  (never,  past,  current,  or  missing),  month  of  blood  draw  (2-month  intervals),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic,  or 
unknown),  and  year  of  diagnosis  (1984-2000  or  2001-2008). 

tModel  further  adjusted  for  time  between  blood  collection  and  cancer  diagnosis  (0  to  <  5  years,  5  to  <  10  years,  or  >  10  years). 

§Model  further  adjusted  for  body  mass  index  (continuous)  and  history  of  diabetes  (yes  or  no). 


engagement  in  stromal  cells.6"8  In  a  study  by  Sherman  and  col¬ 
leagues,8  VDR  acted  as  a  master  transcriptional  regulator  of 
pancreatic  stellate  cells,  which  reside  within  the  tumor  microen¬ 
vironment.  Treatment  of  pancreatic  cancer  in  genetically  engi¬ 
neered  mice  with  a  vitamin  D  analog  led  to  remodeling  of  the 
stroma,  increased  concentration  of  intratumoral  chemotherapy, 
and  prolonged  survival. 

Few  studies  have  investigated  whether  plasma  levels  of  25(OH)D 
predict  patient  survival.  Circulating  25(OH)D  is  a  preferred 
measure  of  vitamin  D  status  as  a  result  of  its  2-week  half-life  and 
its  reflection  of  both  vitamin  D  ingested  in  the  diet  and  synthesized 
in  the  skin.1’  In  a  retrospective  study  of  178  patients  who  were 
treated  at  a  tertiary  cancer  center  and  who  underwent  vitamin  D 
measurement  as  part  of  clinical  care,  serum  25(OH)D  <  20  ng/mL 
at  the  initial  clinical  visit  was  associated  with  poor  prognosis 
(HR,  1.99;  95%  Cl,  1.16  to  3.43)  in  patients  with  advanced 
disease.12  In  contrast,  baseline  serum  25(OH)D  levels  were  not 


Fig  1.  Overall  survival  curves  by  prediagnostic  plasma  25-hydroxyvitamin  D 
[25(OH)D]  among  patients  with  pancreatic  cancer. 


associated  with  overall  survival  in  256  patients  with  advanced 
pancreatic  cancer  who  were  enrolled  in  a  randomized  clinical 
trial  of  gemcitabine  with  or  without  bevacizumab. 1 1  However, 
these  studies  had  relatively  small  sample  sizes,  different  pa¬ 
tient  populations,  and  plasma  25(OH)D  measurements  made 
after  cancer  diagnosis.  Because  inadequate  nutrition  and  lim¬ 
ited  outdoor  activity  as  a  result  of  morbidity  of  recently  diag¬ 
nosed  pancreatic  cancer  can  reduce  25(OH)D  levels,  measured 
levels  likely  do  not  reflect  the  long-term  exposure  to  vitamin 
D  before  acute  illness.  Although  little  data  are  available  related 
to  pancreatic  cancer,  studies  of  several  other  malignancies  have 
suggested  that  single  nucleotide  variants  at  VDR  are  associated  with 
patient  survival,  including  in  smoking-associated  cancers.33"35 

Median  survival  times  were  longer  by  3  to  4  months  in  pa¬ 
tients  with  sufficient  circulating  25(OH)D  compared  with  patients 
with  deficient  levels.  Two  multiagent  chemotherapy  programs  have 
recently  been  adopted  for  treatment  of  patients  with  metastatic 
pancreatic  cancer,  and  have  improved  median  overall  survival  by 
1.8  months  (gemcitabine  plus  nab-paclitaxel)36  and  4.3  months 
(FOLFIRINOX)17  compared  with  single-agent  gemcitabine.  Whether 
addition  of  vitamin  D  or  its  analogs  to  systemic  chemotherapy  can 
further  improve  patient  outcomes  is  currently  being  investigated 
(eg,  ClinicalTrials.gov  NCT02030860). 

Several  strengths  of  this  study  are  notable,  including  large 
sample  size,  inclusion  of  patients  with  all  stages  of  disease  and  from 
all  geographic  regions  of  the  United  States,  and  a  prospective 
cohort  study  design.  An  important  aspect  of  the  prospective  cohort 
design  is  its  ability  to  fully  capture  the  spectrum  of  patients  with 
pancreatic  cancer  in  terms  of  disease  aggressiveness  and  stage  of 
disease,  as  individuals  are  enrolled  before  their  diagnosis  and  are 
not  identified  at  selected  tertiary  care  centers.  Of  note,  survival 
times  and  stage  distribution  for  patients  in  the  five  cohorts  were 
similar  to  121,713  patients  who  were  included  in  the  National 
Cancer  Database,  which  is  thought  to  capture  76%  of  patient  cases 
of  pancreatic  cancer  diagnosed  in  the  United  States  each  year.38  The 
prospective  cohort  design  reduces  bias  that  results  from  reverse 
causation,  as  blood  samples  were  collected  before  nutritional 
deficiencies  and  limited  performance  status  that  commonly  de¬ 
velop  at  the  time  of  pancreatic  cancer  diagnosis.  Furthermore, 
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Table  3. 

HRs  for  Death  Among  Patients  With  Pancreatic  Cancer  by  Quintile  of  Plasma  25(OH)D 

Model 

Quintile  of  Plasma  25(OH)D* 

1 

2 

3 

4 

5 

Median  25(OH)D,  nM 

35.7 

48.2 

58.3 

68.8 

88.1 

Person-months 

1,024 

1,172 

1,371 

1,165 

1,171 

Patients/deaths 

97/92 

99/91 

101/94 

99/95 

97/92 

HR  (95%  Cl)t 

1.0 

0.91  (0.64  to  1.30) 

0.72  (0.51  to  1.01) 

0.74  (0.50  to  1 .09) 

0.68  (0.48  to  0.95) 

HR  (95%  0)4 

1.0 

0.91  (0.64  to  1.31) 

0.77  (0.54  to  1 .09) 

0.76  (0.53  to  1.09) 

0.66  (0.47  to  0.94) 

HR  (95%  CDS 

1.0 

0.96  (0.67  to  1.39) 

0.74  (0.52  to  1 .06) 

0.76  (0.51  to  1.13) 

0.66  (0.46  to  0.95) 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  HR,  hazard  ratio. 

'Cohort-specific  quintile  ranges  of  plasma  25(OH)D:  Health  Professionals  Follow-up  Study  (<  46.9,  46.9  to  <  62.0,  62.0  to  <  72.3,  72.3  to  <  83.9,  a  83.9),  Nurses' 
Health  Study  «  44.3,  44.3  to  <  54. 1 ,  54. 1  to  <  64.7,  64.7  to  <  82.9,  >  82.9),  Physicians'  Health  Study  «  45.0,  45.0  to  <  58.3,  58.3  to  <  68.6,  68.6  to  <  86.9,  >  86.9), 

Women's  Health  Initiative  (<  38.2,  38.2  to  <  49.4,  49.4  to  <  59.7,  59.7  to  <  68.9,  >  68.9),  and  Women's  Health  Study  (<  44.9,  44.9  to  <  56.4,  56.4  to  <  61 .8,  61 .8  to 

<  79.7,  a  79.7). 

tMeta-analysis  of  cohort-specific  HRs  (95%  Cl)  from  Cox  proportional  hazards  regression  models  adjusted  for  age  at  diagnosis,  race  and  ethnicity  (white,  black,  other, 
unknown),  smoking  status  (never,  past,  current,  or  missing),  month  of  blood  draw  (2-month  intervals),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic,  or 
unknown),  and  year  of  diagnosis  (1984-2000  or  2001-2008). 

tModel  further  adjusted  for  time  between  blood  collection  and  cancer  diagnosis  (0  to  <  5  years,  5  to  <  10  years,  or  >  10  years). 

SModel  further  adjusted  for  body  mass  index  (continuous)  and  history  of  diabetes  (yes  or  no). 


exclusion  of  patients  who  were  diagnosed  with  pancreatic  cancer 
within  1  year  of  blood  collection  did  not  materially  alter  our  results. 
Circulating  25(OH)D  was  measured  in  a  single  laboratory  as 
a  single  batch,  with  low  coefficients  of  variance  for  blinded, 
replicate  quality  control  samples,  and  extensive  covariate  data  from 
all  five  cohorts  allowed  for  rigorous  control  of  potential  con¬ 
founding  factors  and  evaluation  of  effect  modification. 

Limitations  of  the  current  study  must  also  be  considered. 
Among  patients  with  pancreatic  cancer,  treatment  programs  likely 
varied,  and  we  could  not  control  for  differences  in  treatment 
because  this  information  was  not  collected  in  our  cohorts.  Nev¬ 
ertheless,  chemotherapy  and  radiation  have  had  only  modest 
impact  on  patient  survival,2  and  treatment  programs  were  unlikely 
to  have  varied  meaningfully  by  baseline  25(OH)D  measured  years 


before  diagnosis.  We  used  overall  mortality  data  in  our  analyses,  as 
opposed  to  pancreatic  cancer-specific  mortality;  however,  pan¬ 
creatic  cancer  is  a  highly  lethal  malignancy,  with  overall  cure  rates 
of  <  5%,  such  that  overall  mortality  is  a  good  surrogate  for  cancer- 
specific  mortality.  Although  circulating  25(OH)D  was  measured  at 
a  single  time  point,  we  have  previously  shown  a  high  correlation  of 
0.70  for  repeated  measures  of  plasma  25(OH)D  within  individuals 
over  time,39  which  suggests  that  a  single  measurement  is  a  rea¬ 
sonable  proxy  for  long-term  levels  of  25(OH)D.  We  cannot  rule  out 
that  our  findings  may  be  influenced,  in  part,  by  residual  con¬ 
founding  or  that  prediagnostic  plasma  25(OH)D  levels  may  mark 
overall  health  status  impacting  patient  survival.  Nonetheless,  we 
included  multiple  possible  confounding  covariates  in  multivariable 
models  without  observing  meaningful  changes  in  risk  estimates. 


Table  4.  HRs  for  Death  Among  Patients  With  Pancreatic  Cancer  by  Plasma  25(OH)D  Level  Stratified  by  Covariates 

Plasma  25(OH)D,  HR  (95%  CD* 


Stratification  Covariate 

No.  of  Patients 

Insufficient 

Relative  Insufficiency 

Sufficient 

P  Interaction 

BMI,  kg/m2 
<  25.0 

216 

1.0 

0.92  (0.63  to  1.35) 

0.86  (0.57  to  1 .30) 

.44 

25.0-29.9 

182 

1.0 

0.67  (0.46  to  0.97) 

0.70  (0.44  to  1.12) 

a  30.0 

95 

1.0 

0.64  (0.35  to  1.17) 

0.37  (0.15  to  0.93) 

Tobacco  use 

Never 

204 

1.0 

0.87  (0.60  to  1.25) 

0.64  (0.41  to  1.01) 

.60 

Past 

221 

1.0 

0.68  (0.48  to  0.96) 

0.68  (0.46  to  1 .00) 

Current 

65 

1.0 

1.12  (0.47  to  2.72) 

1.17  (0.41  to  3.32) 

Cancer  stage 

Localized 

65 

1.0 

0.51  (0.21  to  1.24) 

1.07  (0.39  to  2.96) 

.13 

Locally  advanced 

120 

1.0 

0.91  (0.53  to  1.56) 

0.67  (0.35  to  1 .27) 

Metastatic 

220 

1.0 

0.91  (0.65  to  1.28) 

0.63  (0.41  to  0.96) 

Season  of  blood  draw 

Summer/autumn 

207 

1.0 

0.91  (0.65  to  1.26) 

0.76  (0.48  to  1 .23) 

.66 

Winter/spring 

286 

1.0 

0.71  (0.52  to  0.98) 

0.64  (0.45  to  0.90) 

Time  of  blood  draw  to  diagnosis,  years 
<  5 

165 

1.0 

0.64  (0.43  to  0.95) 

0.58  (0.35  to  0.98) 

.20 

5-10 

183 

1.0 

1.03  (0.70  to  1.50) 

0.75  (0.47  to  1.19) 

>  10 

145 

1.0 

0.88  (0.54  to  1.44) 

1.02  (0.60  to  1.75) 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  BMI,  body  mass  index;  HR,  hazard  ratio. 

'Adjusted  for  age  at  diagnosis,  cohort  (also  adjusts  for  sex),  race  and  ethnicity  (white,  black,  other,  or  unknown),  smoking  status  (never,  past,  current,  or  missing),  month 
of  blood  draw  (2-month  intervals),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic,  or  unknown),  and  year  of  diagnosis  (1984-2000  or  2001-2008). 


www.jco.org  ©  2016  by  American  Society  of  Clinical  Oncology  5 

Information  downloaded  from  jco.ascopubs.org  and  provided  by  at  Harvard  Libraries  on  June  30,  2016  from  128.103.149.52 
Copyright  ©  2016  American  Society  of  Clinical  Oncology.  All  rights  reserved. 


Yuan  et  al 


HPFS  (n  76)  ■ 

HR  (95%  Cl) 

n  89  (n  98  tn  1  88) 

NHS  (n  101)  - ■ - 

0.56(0.30  to  1.01) 

WHI  (n  208)  - ■ - 

0.62  (0.38  to  1.02) 

Meta-analysis 

0.66  (0.49  to  0.90) 

r*°">">gene„y=  ■'»  Q2  Q4  06  Qg 

1.2  1.4 

Fig  2.  Forest  plot  and  meta-analysis  of  hazard  ratios  (HRs)  for  death  among  patients  with  pancreatic  cancer,  comparing  those  with  sufficient  levels  of  25-hydroxyvitamin 
D  with  those  with  insufficient  levels  in  the  Health  Professionals  Follow-up  Study  (HPFS),  Nurses'  Health  Study  (NHS),  Physicians'  Health  Study  (PHS),  Women's  Health 
Initiative  (WHI),  and  Women's  Health  Study  (WHS).  Solid  squares  and  horizontal  lines  correspond  to  the  cohort-specific  multivariable-adjusted  hazard  ratios  and  95%  CIs. 
Area  of  the  solid  square  reflects  the  cohort-specific  weight  (inverse  of  the  variance).  Diamond  represents  the  meta-analysis  multivariable-adjusted  HR  and  95%  Cl.  Vertical 
line  indicates  an  odds  ratio  of  1 .0.  HRs  adjusted  for  age  at  diagnosis,  race  and  ethnicity  (white,  black,  other,  unknown),  smoking  status  (never,  past,  current,  or  missing), 
month  of  blood  draw  (2-month  intervals),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic,  or  unknown),  and  year  of  diagnosis  (1984-2000  or  2001-2008). 


Finally,  our  study  participants  were  predominantly  individuals  of 
European  descent.  African  Americans  have  higher  rates  of  vitamin 
D  deficiency40  and  pancreatic  cancer  mortality.41’44  Additional 
studies  of  vitamin  D  levels  and  pancreatic  cancer  survival  are 
warranted  in  racially  diverse  patient  populations.43 

Higher  prediagnostic  plasma  levels  of  25(OH)D  were  associated 
with  a  statistically  significant  improvement  in  survival  among  patients 
with  pancreatic  cancer  who  were  enrolled  in  five  large,  prospective 
cohorts.  When  considering  these  findings  together  with  previously 
reported  preclinical  data  in  pancreatic  cancer  models,  agonists  of  the 
vitamin  D  receptor  are  a  potentially  attractive  therapeutic  approach 
for  investigation  in  this  highly  lethal  malignancy. 
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Appendix 


Table  A1.  Baseline  Characteristics  of  Patients  With  Pancreatic  Cancer  by  Cohort 


Characteristic 

HPFS 

NHS 

PHS 

WHI 

WHS 

Total 

No.  of  patients 

76 

101 

72 

208 

36 

493 

Plasma  25(OH)D,  nM 

68.5  (24.3) 

65.4  (25.7) 

65.8  (22.1) 

55.2  (19.6) 

60.4  (19.5) 

61.3  (22.7) 

Age  at  blood  draw,  years 

65.6  (7.8) 

60.5  (6.7) 

56.9  (8.6) 

67.0  (7.3) 

58.9  (8.9) 

63.4  (8.5) 

Age  at  diagnosis,  years 

72.7  (8.4) 

70.7  (7.0) 

72.0  (9.4) 

71.9  (7.5) 

63.1  (9.1) 

71.2  (8.3) 

Female  sex 

0  (0) 

101  (100) 

0  (0) 

208  (100) 

36  (100) 

345  (70.0) 

Race/ethnicity 

White 

71  (93.4) 

100  (99.0) 

53  (73.6) 

178  (85.6) 

34  (94.4) 

436  (88.4) 

Black 

1  (1.3) 

1  (1.0) 

1  (1.4) 

14  (6.7) 

1  (2.8) 

18  (3.7) 

Other 

1  (1.3) 

0  (0) 

0  (0) 

14  (6.7) 

0  (0) 

15  (3.0) 

Missing 

3  (3.9) 

0  (0) 

18  (25.0) 

2  (1 .0) 

1  (2.8) 

24  (4.9) 

Body  mass  index,  kg/m2 

25.7  (3.3) 

25.6  (5.0) 

25.7  (2.9) 

27.4  (5.9) 

26.4  (5.3) 

26.4  (5.0) 

Physical  activity,  MET  h/wk 

36.9  (38.5) 

15.4  (17.8) 

12.0  (11.2) 

13.0  (13.2) 

17.8  (21.8) 

17.4  (22.1) 

History  of  diabetes 

3  (3.9) 

5  (5.0) 

3  (4.2) 

15  (7.2) 

3  (8.3) 

29  (5.9) 

Tobacco  use 

Never 

26  (34.2) 

40  (39.6) 

29  (40.3) 

96  (46.2) 

13  (36.1) 

204  (41 .4) 

Past 

42  (55.3) 

41  (40.6) 

30  (41 .7) 

92  (44.2) 

16  (44.4) 

221  (44.8) 

Current 

8  (10.5) 

20  (19.8) 

13  (18.1) 

17  (8.2) 

7  (19.4) 

65  (13.2) 

Missing 

0  (0) 

0  (0) 

0  (0) 

3  (1 .4) 

0  (0) 

3  (0.6) 

Multivitamin  use 

41  (53.9) 

46  (45.5) 

21  (29.2) 

99  (47.6) 

10  (27.8) 

217  (44.0) 

Alcohol  (>  1  drink/d) 

32  (42.1) 

24  (23.8) 

24  (33.3) 

36  (17.3) 

7  (19.4) 

123  (24.9) 

Median  time  blood  draw  to  diagnosis,  years 

6.5 

10.6 

16.6 

5.2 

4.1 

6.7 

Diagnosis  period 

1984-2000 

39  (51.3) 

52  (51.5) 

45  (62.5) 

103  (49.5) 

31  (86.1) 

270  (54.8) 

2001-2008 

37  (48.7) 

49  (48.5) 

27  (37.5) 

105  (50.5) 

5  (13.9) 

223  (45.2) 

Cancer  stage 

Localized 

11  (14.5) 

19  (18.8) 

14  (19.4) 

13  (6.3) 

8  (22.2) 

65  (13.2) 

Locally  advanced 

11  (14.5) 

12  (11.9) 

14  (19.4) 

77  (37.0) 

6  (16.7) 

120  (24.3) 

Metastatic 

33  (43.4) 

44  (43.6) 

32  (44.4) 

91  (43.8) 

20  (55.6) 

220  (44.6) 

Unknown 

21  (27.6) 

26  (25.7) 

12  (16.7) 

27  (13.0) 

2  (5.6) 

88  (17.8) 

Median  survival  time,  months 

All  patients 

5 

5 

6.5 

8 

5 

6 

By  stage 

Localized 

14 

21 

14.5 

25.0 

11.5 

17 

Locally  advanced 

10 

9 

8.5 

12.0 

10.5 

10 

Metastatic 

4 

3 

3.5 

4.0 

4.5 

3 

Unknown 

4 

6 

7.0 

5.0 

3.0 

5 

NOTE.  Continuous  variables  reported  as  mean  (standard  deviation)  and  categorical  variables  reported  as  No.  (%),  unless  otherwise  noted. 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  HPFS,  Health  Professionals  Follow-up  Study;  MET,  metabolic  equivalent;  NHS,  Nurses'  Health  Study;  PHS,  Physicians' 
Health  Study;  WHI,  Women's  Health  Initiative;  WHS,  Women's  Health  Study. 
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Table  A2.  HRs  for  Death  Among  Patients  With  Pancreatic  Cancer  by  Plasma  25(OH)D  Level,  Excluding  Patients  Diagnosed  Within  1  Year  of  Blood  Draw 

Model 

Plasma  25(OH)D 

^trend 

Insufficient 

Relative  Insufficiency 

Sufficient 

Range 

ng/mL 

<  20 

20  to  <  30 

o 

CO 

Al 

nM 

<  50 

50  to  <  75 

-  75 

Median  survival,  months 

5 

7 

8 

Person-months 

1,662 

2,541 

1,325 

Patients/deaths 

159/148 

201/1 89 

113/108 

HR  (95%  Cl)t 

1.0 

0.84  (0.53  to  1.33) 

0.62  (0.45  to  0.85) 

.02 

HR  (95%  0)4 

1.0 

0.87  (0.55  to  1.36) 

0.63  (0.45  to  0.87) 

.03 

HR  (95%  CDS 

1.0 

0.75  (0.41  to  1.36) 

0.59  (0.41  to  0.84) 

.04 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  HR,  hazard  ratio. 

*Test  for  trend  calculated  by  entering  log-transformed  25(OH)D  as  a  continuous  variable  in  Cox  proportional  hazards  regression  models. 

t  Meta-analysis  of  cohort-specific  HRs  (95%  Cl)  from  Cox  proportional  hazards  regression  models  adjusted  forage  at  diagnosis,  race  and  ethnicity  (white,  black,  other,  or 
unknown),  smoking  status  (never,  past,  current,  or  missing),  month  of  blood  draw  (2-month  intervals),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic,  or 
unknown),  and  year  of  diagnosis  (1984-2000  or  2001-2008). 

4Model  further  adjusted  for  time  between  blood  collection  and  cancer  diagnosis  (0  to  <  5  years,  5  to  <  10  years,  or  >  10  years). 

§  Model  further  adjusted  for  body  mass  index  (continuous)  and  history  of  diabetes  (yes  or  no). 


Table  A3.  SNPs  in  VDR  Ranked  by  Association  P  Value  With  Overall  Survival  Among  Patients  With  Pancreatic  Cancer 

Additive  Model 


SNP* * 

No.  of  Patients 

Minor  Allele  Frequency,  % 

HR  (95%  Cl) 

P 

rs7299460 

391 

0.30 

0.80  (0.68  to  0.95) 

.01 

rsl  1 568820 

388 

0.22 

0.79  (0.65  to  0.96) 

.02 

rs4334089 

388 

0.23 

0.81  (0.68  to  0.98) 

.03 

rs886441 

392 

0.18 

0.82  (0.67  to  0.99) 

.04 

rs2853564 

396 

0.38 

1.16  (0.99  to  1.35) 

.07 

rs22391 79 

392 

0.44 

1.14  (0.99  to  1.33) 

.08 

rs22391 82 

392 

0.47 

0.89  (0.77  to  1 .03) 

.12 

rs7295021 

395 

0.20 

0.86  (0.70  to  1 .05) 

.13 

rsl  1168293 

392 

0.34 

1.13  (0.96  to  1.32) 

.15 

rs22391 86 

395 

0.20 

0.89  (0.74  to  1 .06) 

.20 

rs71 32324 

391 

0.35 

1.10  (0.94  to  1.29) 

.24 

rs7963776 

388 

0.46 

0.91  (0.79  to  1 .06) 

.24 

rs2254210 

396 

0.39 

1.09  (0.94  to  1.26) 

.28 

rs21 89480 

390 

0.34 

0.92  (0.79  to  1 .08) 

.30 

rs2283342 

392 

0.15 

0.90  (0.72  to  1.11) 

.31 

rs4760648 

391 

0.43 

0.92  (0.78  to  1 .08) 

.32 

rsl  0875693 

389 

0.33 

1.08  (0.92  to  1.27) 

.35 

rs2238136 

393 

0.27 

0.92  (0.78  to  1 .09) 

.35 

rsl  1 574077 

391 

0.05 

1.18  (0.83  to  1.66) 

.36 

rsl  859281 

395 

0.07 

1.13  (0.83  to  1.53) 

.44 

rsl  544410 

359 

0.37 

1.05  (0.90  to  1.23) 

.53 

rs7310552 

393 

0.40 

1.04  (0.88  to  1.22) 

.64 

rsl  0747524 

329 

0.49 

0.97  (0.81  to  1.16) 

.71 

rs21 07301 

394 

0.28 

0.97  (0.83  to  1.14) 

.71 

rsl  1168275 

392 

0.24 

0.97  (0.81  to  1.16) 

.73 

rs2544037 

394 

0.42 

1.02  (0.87  to  1.18) 

.84 

rsl  1 574032 

297 

0.09 

1.02  (0.75  to  1.40) 

.89 

rsl  2721 364 

391 

0.14 

1.01  (0.81  to  1.26) 

.93 

rsl  1574143 

391 

0.11 

1.01  (0.80  to  1.27) 

.96 

rs3819545 

393 

0.39 

1.00  (0.86  to  1.17) 

.96 

Abbreviations:  HR,  hazard  ratio;  SNP,  single  nucleotide  polymorphism;  VDR,  vitamin  D  receptor. 

*Effect  of  each  minor  allele  of  the  SNP  on  survival  in  Cox  proportional  hazards  regression  models  adjusted  for  age  at  diagnosis  (years,  continuous),  cohort  (also  adjusts 
for  sex),  race  and  ethnicity  (white,  black,  other,  or  unknown),  smoking  status  (never,  past,  current,  or  missing),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic, 
or  unknown),  and  year  of  diagnosis  (1984-2000  or  2001-2008). 
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Table  A4.  HRs  (95%  Cl)  for  Death  Among  Patients  With  Pancreatic  Cancer  by  Plasma  25(OH)D  and  Stratified  by  rs7299460  Genotypes  in  1 /DR 

Plasma  25(OH)D,  HR  (95%  Cl) 


SNP 

Genotype 

No. 

Insufficient 

Relative  Insufficiency 

Sufficient 

rs7299460 

GG 

194 

1.0 

0.97  (0.69  to  1.38) 

0.78  (0.51  to  1.18) 

AG 

162 

0.84  (0.57  to  1.25) 

0.67  (0.46  to  0.98) 

0.85  (0.56  to  1 .30) 

AA 

35 

0.71  (0.31  to  1.59) 

0.62  (0.33  to  1.16) 

0.45  (0.22  to  0.93) 

NOTE.  Adjusted  for  age  at  diagnosis  (years,  continuous),  cohort  (also  adjusts  for  sex),  race  and  ethnicity  (white,  black,  other,  or  unknown),  smoking  status  (never,  past, 
current,  or  missing),  month  of  blood  draw  (2-month  intervals),  stage  at  diagnosis  (localized,  locally  advanced,  metastatic,  or  unknown),  and  year  of  diagnosis  (1 984-2000 
or  2001-2008). 

Abbreviations:  25(OH)D,  25-hydroxyvitamin  D;  SNP,  single  nucleotide  polymorphism;  VDR,  vitamin  D  receptor. 
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ABSTRACT 


Leptin  is  an  adipokine  involved  in  regulating  energy  balance  which  has  been  identified  as  a 
potential  biological  link  in  development  of  obesity-associated  cancers,  such  as  pancreatic  cancer. 
In  this  prospective,  nested  case-control  study  of  470  cases  and  1094  controls  from  five  U.S. 
cohorts,  we  used  conditional  logistic  regression  to  evaluate  pancreatic  cancer  risk  by 
prediagnostic  plasma  leptin,  adjusting  for  race/ethnicity,  diabetes,  body-mass  index,  physical 
activity,  plasma  C-peptide,  adiponectin  and  2 5 -hydroxy vitamin  D.  Due  to  known  differences  in 
leptin  levels  by  gender,  analyses  were  conducted  separately  for  men  and  women.  We  also 
evaluated  associations  between  32  tagging  single  nucleotide  polymorphisms  (SNPs)  in  the  leptin 
receptor  ( LEPR )  gene  and  pancreatic  cancer  risk.  Leptin  levels  were  higher  in  female  versus 
male  control  participants  (median,  20.8  vs.  6.7pg/mL;  PO.OOO  1 ).  Among  men,  plasma  leptin 
was  positively  associated  with  pancreatic  cancer  risk,  and  those  in  the  top  quintile  had  a 
multivariable-adjusted  odds  ratio  (OR)  of  3.02  (95%  Cl,  1.27-7.16;  /Vtk1=0.02)  compared  to  men 
in  the  bottom  quintile.  Among  women,  circulating  leptin  was  not  associated  with  pancreatic 
cancer  risk  (7Vnd=0.2 1 ).  Results  were  similar  across  cohorts  (Pheterogeneity=0.88  for  two  male 
cohorts  and  0.35  for  three  female  cohorts).  In  genetic  analyses,  rsl0493380  in  LEPR  was 
associated  with  increased  pancreatic  cancer  risk  among  women,  with  an  OR  per  minor  allele  of 
1.54  (95%  Cl,  1.18-2.02;  multiple  hypothesis-corrected  P=0.03).  No  SNPs  were  significantly 
associated  with  risk  in  men.  In  conclusion,  higher  prediagnostic  levels  of  plasma  leptin  were 
associated  with  an  elevated  risk  of  pancreatic  cancer  among  men,  but  not  among  women. 
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INTRODUCTION 


Pancreatic  cancer  is  the  third-leading  cause  of  cancer-related  death  in  the  United  States 
(1).  Obesity  is  associated  with  increased  pancreatic  cancer  risk  (2,3),  but  the  underlying 
mechanisms  are  poorly  understood.  Leptin  was  initially  discovered  in  studies  of  obese  mice, 
which  were  noted  to  have  increased  caloric  intake.  Subsequent  studies  demonstrated  that  the 
underlying  cause  of  the  obese  phenotype  in  these  mice  was  a  truncating  mutation  in  the  gene 
encoding  for  leptin  (4).  Leptin  is  a  hormone  synthesized  almost  exclusively  in  adipocytes,  and 
plasma  levels  of  leptin  in  humans  are  proportional  to  total  body  adipose  tissue  (5).  When 
circulating  levels  of  leptin  are  increased,  downstream  signaling  is  activated  through  the  ObR 
(HUGO  Gene  Nomenclature  Committee  [HGNC]  Symbol  LEPR)  transmembrane  receptor  on 
target  cells  (6).  The  leptin  receptor  is  present  on  cells  within  the  hypothalamus,  which  is  the 
mechanism  by  which  leptin  is  thought  to  regulate  caloric  intake  (7).  Additionally,  leptin  receptor 
expression  is  distributed  widely  throughout  the  body,  including  on  cells  within  the  pancreas  (8). 
Leptin  receptors  are  also  present  on  tumor  cells  and  activation  of  these  receptors  increases  cancer 
cell  proliferation  and  reduces  rates  of  cancer  cell  apoptosis  (9). 

Given  its  essential  role  in  regulating  energy  balance,  leptin  may  be  an  important 
biological  link  in  the  development  of  obesity-associated  malignancies,  including  pancreatic 
cancer.  To  further  investigate  the  role  of  leptin  in  pancreatic  carcinogenesis,  we  examined  the 
association  between  prediagnostic  plasma  leptin  and  subsequent  pancreatic  cancer  risk  in  five 
prospective  cohorts  with  up  to  26  years  of  follow-up  since  blood  collection.  We  additionally 
evaluated  the  association  between  polymorphisms  in  LEPR  with  risk  of  pancreatic  cancer. 
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METHODS 


Study  Participants 

We  pooled  blood  samples  and  data  from  five  U.S.  prospective  cohorts.  The  Health 
Professionals  Follow-up  Study  (HPFS)  enrolled  51,529  male  health  professionals  aged  40-75  in 
1986.  The  Nurses’  Health  Study  (NHS)  enrolled  121,700  female  nurses  aged  30-55  in  1976.  The 
Physicians’  Health  Study  I  (PHS)  was  a  randomized  clinical  trial  of  aspirin  and  P-carotene  that 
enrolled  22,071  healthy  male  physicians  aged  40-84  in  1982.  The  aspirin  component  of  the  trial 
ended  in  1988,  while  the  P-carotene  component  ended  in  1995,  and  PHS  I  participants  continue 
follow-up  as  an  observational  cohort.  The  Women’s  Health  Initiative  (WHI)-Observational 
Study  enrolled  93,676  postmenopausal  women  aged  50-79  between  1994  and  1998.  The 
Women’s  Health  Study  (WHS)  was  a  randomized  clinical  trial  of  low-dose  aspirin  and  vitamin  E 
that  enrolled  39,876  healthy  female  health  professionals  aged  >45  in  1992.  The  trial  was 
completed  in  2004  and  WHS  participants  continue  to  be  followed  as  an  observational  cohort. 

Individual  characteristics  and  lifestyle  factors  were  obtained  from  baseline  questionnaires 
at  enrollment  in  PHS,  WHI,  and  WHS  and  from  the  questionnaires  preceding  blood  draw  in 
HPFS  and  NHS.  Details  of  these  cohorts  have  been  described  previously  (10-14).  The  current 
study  was  approved  by  the  Human  Research  Committee  at  the  Brigham  and  Women’s  Hospital, 
Boston,  MA,  and  participants  provided  informed  consent. 

Blood  Collection  and  Plasma  Assays 

Blood  samples  were  collected  from  18,225  men  in  HPFS  from  1993-1995,  32,826 
women  in  NHS  from  1989-1990,  14,916  men  in  PHS  from  1982-1984,  93,676  women  in  WHI 
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from  1994-1998,  and  28,345  women  in  WHS  from  1992-1995.  Details  on  blood  draw, 
transportation,  and  storage  have  been  described  previously  (12,14-16). 

Plasma  leptin  was  assayed  in  the  laboratory  of  Dr.  Nader  Rifai  (Children’s  Hospital, 
Boston,  MA),  using  reagents  from  R&D  Systems  (Minneapolis,  MN).  Measurement  of  plasma 
adiponectin  and  C-peptide  were  described  previously  (17)  .  All  samples  for  leptin  and 
adiponectin  were  handled  identically  in  a  single  batch  and  C-peptide  was  handled  in  two  batches. 
Laboratory  personnel  were  blinded  to  case  or  control  status.  The  mean  intra-assay  coefficients  of 
variance  for  quality  control  samples  were  <10%  for  each  biomarker. 

Pancreatic  cancer  cases  and  matched  controls 

We  included  cases  of  pancreatic  adenocarcinoma  diagnosed  through  2008  with 
prediagnostic  blood  and  no  prior  history  of  cancer,  except  non-melanoma  skin  cancer.  Incident 
cases  were  identified  by  self-report  or  during  follow-up  of  a  participant’s  death.  Deaths  were 
ascertained  from  next-of-kin  or  the  U.S.  postal  service  and  by  searching  the  National  Death 
Index.  Medical  records  of  the  cases  were  requested  and  reviewed  by  study  physicians  blinded  to 
exposure  data. 

Eligible  controls  were  cohort  participants  who  provided  blood  and  were  alive  and  free  of 
cancer  at  the  date  of  the  case’s  diagnosis.  We  randomly  selected  2-3  controls  for  each  case, 
matching  on  year  of  birth  (±  5  years),  prospective  cohort  (which  concurrently  matched  on  sex), 
smoking  status  (never,  past,  current),  fasting  status  (fasting,  non-fasting),  and  month/year  of 
blood  draw. 

For  the  present  analysis,  488  pancreatic  cancer  cases  and  1132  matched  controls  with 
plasma  were  available.  Due  to  concern  regarding  the  possible  influence  of  subclinical 
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malignancy,  we  excluded  pancreatic  cancer  cases  diagnosed  within  1  year  of  blood  draw  (n=19) 
and  their  matched  controls  (n=38),  resulting  in  a  total  of  470  cases  and  1094  controls 
(Supplementary  Table  1).  Of  these  470  cases,  465  (99%)  were  confirmed  by  review  of  medical 
records,  tumor  registry  data,  or  death  certificates. 

Single  Nucleotide  Polymorphisms  (SNP)  selection  and  Genotyping 

A  total  of  39  SNPs  in  the  LEPR  gene  +/-  20kb  were  selected  with  the  tagger  algorithm  in 
Haploview,  using  r=0.8  and  minor  allele  frequency  (MAF)  >5%  among  Whites  from  the 
HapMap  Project  database.  Five  SNPs  associated  with:  serum  amyloid  A  (rsl275319,  (18)), 
soluble  Ob-R  (rs2767485,  (19)),  and  CRP:  (rs4420065  (20),  rs6700896  (21)  and  rsl892534  (22) 
were  forced  in.  From  412  cases  (Supplementary  Table  1),  DNA  was  extracted  from  huffy  coat 
using  QIAGEN  QIAmp  and  whole  genome  amplified  using  GE  Healthcare  GenomiPhi. 
Genotyping  was  performed  at  Partners  Healthcare  Center  for  Personalized  Genetic  Medicine 
using  a  custom-designed  Illumina  Golden  Gate  genotyping  assay.  Seven  tagging  SNPs  were  not 
supported  by  the  platfonn.  One  SNP  (rs913199)  deviated  from  Hardy- Weinberg  Equilibrium  at 
P=0.008.  Replicate  samples  tested  for  quality  control  (N=44  groups)  had  a  mean  genotype 
concordance  of  98.2%. 

Statistical  Analysis 

Median  leptin  levels  among  cases  and  controls  were  compared  using  the  Wilcoxon  rank- 
sum  test.  Since  men  and  women  have  different  distributions  of  leptin  levels  (23),  we  performed 
separate  analyses  by  gender  using  pooled  gender-specific  quintiles  from  controls.  To  compute 
odds  ratios  (ORs)  and  95%  confidence  intervals  (CIs),  we  used  conditional  logistic  regression.  In 
multivariate  models,  we  adjusted  for  potential  confounding  factors,  including  race  (White,  Black, 
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other),  multivitamin  use  (yes,  no),  diabetes  (yes,  no),  body  mass  index  (BMI,  kg/nr),  physical 
activity  (MET-hr/wk),  plasma  C-peptide  (continuous),  plasma  adiponectin  (quartiles,  as  previous 
analysis  demonstrated  non-linear  association  of  adiponectin  and  pancreatic  cancer  risk  (17)),  and 
plasma  25(OH)D  (continuous),  /-’-trends  were  calculated  by  the  Wald  test  of  a  score  variable  that 
contained  median  values  of  quintiles.  We  also  conducted  a  meta-analysis  of  cohort-level  data 
among  men  and  women.  We  calculated  ORs  for  each  cohort  and  then  pooled  the  ORs  to  compute 
a  summary  OR  by  gender  using  the  random  effects  model  (24).  Heterogeneity  across  studies  was 
tested  using  the  Q  statistic  (24).  To  evaluate  whether  the  association  between  leptin  and 
pancreatic  cancer  risk  was  log-linear,  we  compared  the  model  fit  including  linear  and  cubic 
spline  terms  to  the  model  fit  with  only  the  linear  tenn  using  the  likelihood  ratio  test  (25).  We 
conducted  subgroup  analyses  using  unconditional  logistic  regression  adjusted  for  the  matching 
factors  and  covariates.  Tests  for  interaction  were  performed  by  the  Wald  test  of  cross-product 
tenns.  We  conducted  sensitivity  analyses  excluding  diabetics  or  cases  diagnosed  within  2  or  4 
years  from  blood  draw. 

We  examined  the  association  between  LEPR  SNPs  and  pancreatic  cancer  risk  by 
modeling  each  genotype  as  number  of  copies  of  the  minor  allele  (additive  model)  using 
conditional  logistic  regression.  We  used  R  software  (version  R.3.2.2.)  to  calculate  the  corrected 
P- value  by  taking  into  account  the  total  number  of  comparisons,  as  well  as  correlations  between 
32  SNPs(26).  We  used  HaploReg  v4.1  to  explore  the  noncoding  functional  characteristics  of 
identified  and  highly  correlated  SNPs  (r  >0.6  in  1000G  CEU  data).  Statistical  analyses  were 
performed  with  SAS  9.1  (SAS  Institute,  Cary,  North  Carolina),  and  all  P  values  are  two  sided. 


RESULTS 
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The  median  time  between  blood  collection  and  cancer  diagnosis  was  7. 1  years  among 
cases.  Among  controls,  median  plasma  leptin  was  20.8  pg/mL  for  women  and  6.7  pg/mL  for 
men.  Leptin  levels  were  comparable  across  studies  for  men  (HPFS  and  PHS)  and  for  women 
(NHS,  WHI,  and  WHS)  (Supplementary  Table  2).  Individuals  with  higher  leptin  levels  had 
higher  BMI,  plasma  C-peptide,  and  prevalence  of  diabetes  (Table  1).  After  adjusting  for  age, 
fasting  status  and  cohort,  Spearman  correlation  coefficients  for  plasma  leptin  and  BMI  were  0.50 
(P<0.0001)  among  men  and  0.73  (P<0.000 1 )  among  women  (Figure  1)  (Supplementary  table  3), 
similar  to  those  reported  in  other  populations  (27,28). 

We  observed  a  positive  association  between  plasma  leptin  and  pancreatic  cancer  risk 
among  men,  but  not  among  women  (Pheterogeneity=0.02;  Table  2).  In  the  base  model  conditioned  on 
matching  factors,  compared  to  the  bottom  quintile,  men  in  the  top  quintile  had  an  OR  of  2.77 
(1.37-5.61),  (Ptrend=0.01;  Table  2).  In  comparison,  women  in  the  top  quintile  had  an  OR  of  1.27 
(0.84-1.91)  (Ptrend=0.64;  Table  2).  Further  adjustment  for  race,  multivitamin  use,  plasma 
25(OH)D,  history  of  diabetes,  BMI,  physical  activity,  plasma  C-peptide,  and  plasma  adiponectin 
yielded  similar  results  (Table  2).  Similar  associations  were  observed  in  sensitivity  analyses  when 
we  excluded  cases  with  diabetes  or  cases  diagnosed  within  2  or  4  years  of  blood  collection  and 
their  matched  controls  (Supplementary  Table  4). 

Spline  curves  were  consistent  with  log-linear  associations  (Pnoniinear=0.81  for  men; 

P nOniinear=0- 14  for  women).  Therefore,  in  subsequent  meta-analyses  and  subgroup  analyses,  we 
modeled  leptin  as  a  continuous  variable.  The  multivariate  ORs  for  an  increment  of  5  ng/mL  in 
plasma  leptin  were  1.25  (1.02-1.54)  for  men  and  0.98  (0.93-1.03)  for  women  (Supplementary 
Table  5).  ORs  were  similar  within  the  two  male  and  the  three  female  cohorts  (Figure  2; 
Rheterogeneity=0.88  for  HPFS  and  PHS,  and  P heterogeneity=0 . 3 5  for  NHS,  WHI,  and  WHS).  In 
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stratified  analyses,  no  statistically  significant  effect  modification  was  observed  (Supplementary 
Table  5). 

Several  SNPs  at  the  LEPR  gene  were  associated  with  pancreatic  cancer  risk  among 
women  to  PA).  05  in  an  additive  genetic  model  (Table  3,  Supplementary  Figure  1).  After 
adjusting  for  multiple  comparisons,  rsl0493380  located  intronic  to  LEPR  remained  statistically 
significantly  associated  with  increased  risk  of  pancreatic  cancer  (OR  per  minor  allele=1.54;  95% 
CI=1.18-2.02,  multiple  hypothesis-corrected  P=0.03;  Table  3).  This  association  was  consistent 
across  the  three  female  cohorts  (Pheterogeneity=0.28).  The  association  for  rsl0493380  was  not 
statistically  significant  among  men  (OR=1.19,  95%  01=0.79-1.78,  multiple  hypothesis-corrected 
P=1.00;  Supplementary  Table  6).  Analysis  of  rsl0493380  and  highly  correlated  SNPs  using 
HaploReg  identified  multiple  transcription  factor  binding  altered  by  these  SNPs  (Supplementary 
Table  7).  Furthermore,  in  a  blood  eQTL  database  (29),  rsl0493380  (index  SNP)  and  rs3790429 
(in  high  LD  with  rsl0493380,  r2=0.90  in  1000G  EUR)  were  found  to  have  cis  eQTL  effects  on 
LEPR  gene  expression.  No  statistically  significant  association  was  identified  between  SNPs  at 
the  LEPR  gene  and  plasma  leptin  levels  among  controls  (Supplementary  Table  8). 
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DISCUSSION 


As  seen  in  prior  studies  (23),  plasma  leptin  levels  were  higher  in  our  female  compared 
with  male  control  subjects,  with  a  median  level  approximately  three  times  higher  in  women 
versus  men.  Interestingly,  higher  prediagnostic  plasma  leptin  was  associated  with  an  increased 
risk  of  pancreatic  cancer  in  men,  with  no  increase  in  risk  was  observed  in  women.  This  positive 
association  among  men  was  independent  of  other  known  risk  factors  for  pancreatic  cancer, 
including  characteristics  and  plasma  markers  associated  with  obesity  and  insulin  resistance. 
Furthermore,  the  association  was  highly  consistent  across  two  cohorts  with  male  participants 
(positive  association  in  HPFS  and  PHS)  and  three  cohorts  with  female  participants  (no 
association  in  NHS,  WHI,  and  WHS).  In  contrast,  single  nucleotide  variants  at  the  leptin  receptor 
(LEPR)  gene  were  associated  with  pancreatic  cancer  risk  only  in  women,  but  not  in  men.  This 
association  was  consistent  across  the  three  female  cohorts.  Notably,  the  most  statistically 
significant  SNP  (rsl0493380)  and  highly  correlated  variants  at  LEPR  may  alter  leptin  receptor 
gene  expression  based  on  bioinformatic  analyses.  In  aggregate,  these  data  support  the  importance 
of  adipokines  and  adipokine  signaling  in  pancreatic  cancer  risk  in  men  and  women,  even  though 
baseline  circulating  leptin  levels  differ  greatly  by  gender. 

A  previous  nested  case-control  study  evaluated  the  association  of  prediagnostic  plasma 
leptin  with  risk  of  pancreatic  cancer  (30).  In  this  pooled  analysis  of  three  cohorts,  plasma  leptin 
was  not  associated  with  risk  of  pancreatic  cancer  during  the  first  10  years  of  follow-up,  using 
gender-specific  quintiles  of  plasma  leptin.  However,  a  statistically  significant  positive 
association  was  observed  among  men  (OR,  2.94;  95%  Cl,  1.64-6.46;  Ptrend=0.001;  comparing 
extreme  quintiles)  with  longer  follow-up  time  (>10  years),  while  the  association  among  women 
for  this  duration  could  not  be  evaluated  due  to  small  sample  size.  Similarly,  in  our  cohorts,  a 
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positive  association  of  plasma  leptin  with  pancreatic  cancer  risk  was  observed  in  men.  However, 
in  stratified  analyses  by  time  between  blood  collection  and  cancer  diagnosis,  statistically 
significant  effect  modification  was  not  seen.  This  difference  may  be  partly  due  to  our  a  priori 
exclusion  of  cases  with  blood  collected  within  12  months  of  cancer  diagnosis,  the  time  period 
during  which  cancer-associated  weight  loss  most  commonly  occurs  (31).  The  stratified  analyses 
in  men  were  also  limited  by  smaller  sample  sizes  within  strata.  Three  previous  cross-sectional 
studies  observed  that  plasma  leptin  levels  were  lower  in  pancreatic  cancer  patients  compared  to 
controls  (32-34).  However,  in  these  retrospective  studies,  hypoleptinemia  may  have  been  due  to 
the  weight  loss  that  is  commonly  experienced  by  patients  with  pancreatic  cancer  (5).  Therefore, 
it  is  difficult  to  determine  whether  the  observed  low  leptin  levels  contributed  to  pancreatic 
carcinogenesis  or  were  a  consequence  of  the  cancer. 

Several  lines  of  evidence  support  a  biological  link  between  leptin  and  pancreatic 
carcinogenesis.  Leptin  plays  a  central  role  in  the  regulation  of  insulin  sensitivity  (35),  and  studies 
have  demonstrated  associations  between  hyperglycemia,  insulin  resistance,  and  future  risk  of 
pancreatic  cancer  (17,36-38).  Therefore,  one  mechanism  by  which  leptin  may  influence 
pancreatic  cancer  risk  is  through  its  modulation  of  insulin  sensitivity  (39).  Leptin  is  synthesized 
by  adipose  tissue,  and  its  concentration  correlates  with  total  body  fat.  Leptin  may  therefore  act  as 
a  better  marker  for  the  relevant  states  of  adiposity  than  BMI,  which  does  not  discriminate 
between  fat  and  muscle  mass  (40).  Target  tissue  effects  of  circulating  leptin  are  not  solely 
mediated  centrally  on  cells  within  the  hypothalamus,  as  leptin  receptors  are  widely  distributed  in 
the  body  (7).  Leptin  receptors  have  been  identified  on  the  surface  of  tumor  cells,  including 
pluripotent  cells  thought  to  function  as  tumor  initiating  cells  (41).  Therefore,  leptin  signaling 
may  have  an  important  direct  role  in  promoting  tumor  initiation  and  growth,  independent  of  its 
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role  in  insulin  sensitivity  (42).  As  a  consequence  of  this  potential  direct  effect,  inhibitors  of  the 
leptin  receptor  are  being  explored  as  novel  therapeutics  to  inhibit  tumor  growth  in  patients  (43). 

We  observed  a  positive  association  between  plasma  leptin  and  pancreatic  cancer  risk  only 
in  men.  Interestingly,  leptin  levels  are  considerably  lower  in  men  than  in  women,  even  for  the 
same  age  and  body-mass  index  (44).  Although  the  underlying  reasons  for  this  difference  in 
circulating  leptin  are  unclear,  sex-differences  in  reproductive  hormones  and  body  fat  distribution 
have  been  proposed  as  possible  etiologies.  Particularly,  women  tend  to  have  higher  total  and 
subcutaneous  fat,  while  men  have  a  greater  percentage  of  visceral  fat,  which  may  influence 
circulating  levels  of  leptin  (40).  Given  the  substantially  lower  leptin  levels  and  differing 
metabolic  environment  in  men,  the  actions  of  leptin  may  be  sex-specific,  with  implications  for 
disease  development.  Alternatively,  physiologic  differences  related  to  circulating  leptin  may  be 
detectable  only  when  compared  in  the  lower  ranges  of  circulating  leptin,  which  are  seen 
predominantly  in  men.  Notably,  a  number  of  prospective  cohort  studies  have  examined  the 
association  of  prediagnostic  plasma  leptin  and  the  subsequent  diabetes  risk.  Similar  to  the  current 
study  of  pancreatic  cancer,  most  of  these  studies  demonstrated  positive  associations  of  leptin  and 
diabetes  in  men,  but  not  in  women  (45).  Of  note,  crosstalk  between  estrogen  and  leptin  signaling 
has  been  shown  previously  (46).  Other  than  modulating  synthesis  of  leptin  (47),  and  leptin 
receptor  (48),  estrogen  receptor  alpha  enhances  leptin-induced  activation  of  downstream 
signaling  pathways,  including  the  JAK/STAT  pathway  (46).  Alternatively,  the  different 
association  between  leptin  levels  and  pancreatic  cancer  risk  between  genders  in  our  study  could 
be  due  to  chance,  and  these  findings  should  be  confirmed  in  future  studies. 

The  current  study  has  several  strengths.  The  prospective  design  and  exclusion  of  cases 
diagnosed  within  1  year  of  blood  collection  reduced  the  potential  impact  of  reverse  causation  on 
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our  results.  Furthermore,  similar  associations  were  observed  when  we  excluded  cases  diagnosed 
within  2  or  4  years  of  blood  collection.  Leptin  was  measured  in  a  single  laboratory  as  a  single 
batch,  with  low  coefficients  of  variance  for  quality  control  samples.  In  our  analyses,  we  adjusted 
for  BMI,  physical  activity,  and  other  biomarkers  related  to  insulin  resistance,  including  C- 
peptide  and  adiponectin,  to  rigorously  control  for  confounding.  We  evaluated  not  only 
circulating  leptin  levels,  but  also  genetic  variants  in  the  leptin  receptor,  which  may  affect  signal 
transduction  after  ligand  binding.  Additional  strengths  included  a  large  sample  size,  long  follow¬ 
up  period,  and  inclusion  of  men  and  women. 

Plasma  leptin  was  measured  at  only  one  point  in  time,  so  leptin  levels  may  not  fully 
reflect  long-tenn  plasma  concentrations.  However,  leptin  levels  are  relatively  stable  over  time  in 
healthy  subjects;  repeated  plasma  leptin  measurements  1-year  apart  demonstrated  a  high  intra¬ 
class  correlation  coefficient  of  0.74  (21).  We  cannot  rule  out  that  residual  confounding  by 
adiposity  not  captured  by  BMI  may  be  present;  however,  adjustment  for  plasma  adiponectin, 
25(OH)D  and  C-peptide  as  markers  of  adiposity  and  insulin  resistance  did  not  materially  alter 
our  results.  We  identified  an  association  of  rsl0493380  with  pancreatic  cancer  risk  in  women. 
This  SNP  was  not  identified  as  genome-wide  significant  in  two  recent  pancreatic  cancer  genome¬ 
wide  association  studies  (49,50).  However,  these  studies  were  required  to  meet  a  stringent 
multiple-hypothesis  testing  threshold  for  statistical  significance  due  to  testing  of  >500,000  SNPs, 
included  approximately  60%  men,  and  included  a  majority  of  patients  from  tertiary  center  cases- 
control  studies.  Finally,  our  study  population  consisted  primarily  of  White  participants  and 
further  studies  are  required  to  evaluate  circulating  leptin  in  participants  of  different 
race/ethnicity. 
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In  conclusion,  we  identified  a  positive  association  between  prediagnostic  circulating 
leptin  levels  and  pancreatic  cancer  risk  in  men,  independent  of  obesity  and  other  risk  factors  of 
pancreatic  cancer.  Although  this  association  was  not  observed  in  women,  single  nucleotide 
variants  in  the  leptin  receptor  were  associated  with  pancreatic  cancer  risk  in  women,  and 
bioinformatic  analyses  suggested  differences  in  leptin  receptor  expression  with  these  variants. 
Our  data  provide  additional  evidence  for  a  biological  link  between  obesity,  insulin  resistance, 
and  pancreatic  cancer  risk,  specifically  focusing  attention  on  adipokines  and  adipokine  signaling 
in  pancreatic  cancer  development. 
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Table  1.  Age-  and  study-standardized  baseline  characteristics  according  to  plasma  leptin  levels 
among  controls 


Quintiles  of  plasma  leptin 


Characteristic* 

1 

2 

3 

4 

5 

Men 

Plasma  leptin  level,  ng/mL 

<3.0 

3. 1-4.8 

4.9-7. 8 

7.9-11.7 

>11.8 

No.  of  controls 

69 

71 

71 

72 

70 

Age  at  blood  draw,  years 

60.7(10.1) 

60.1  (8.9) 

60.2(10.0) 

60.8  (8.9) 

61.0(9.1) 

Race,  % 

White 

94.0 

93.2 

85.5 

93.7 

91.3 

Black 

0 

1.0 

0 

0 

0 

Other 

6.0 

5.9 

14.5 

6.3 

8.7 

Body-mass  index,  kg/m2 

23.4  (2.6) 

24.7  (2.1) 

24.8  (2.3) 

26.0  (2.5) 

27.6  (3.6) 

Physical  activity,  MET- 
hr/week 

34.6  (34.9) 

29.1  (31.6) 

30.6  (48.6) 

27.2  (34.6) 

18.1  (22.6) 

Cigarette  smoking,  % 

Never 

37.0 

36.9 

37.3 

40.9 

33.6 

Past 

52.9 

50.4 

40.7 

41.1 

47.8 

Current 

10.2 

12.7 

22.0 

18.1 

18.6 

Missing 

0 

0 

0 

0 

0 

History  of  diabetes  mellitus, 

% 

1.0 

2.5 

2.2 

1.6 

7.3 

Regular  multivitamin  use, 

% 

37.8 

37.1 

36.0 

30.2 

30.2 

Plasma  25(OH)D,  n mo  1/1. 

76.9  (28.9) 

74.5  (32.0) 

72.3  (20.7) 

70.2  (20.9) 

66.3  (19.5) 

Plasma  C-peptide  levels, 
ng/mL 

1.8  (1.1) 

2.2  (1.4) 

2.5  (1.6) 

2.8  (3.4) 

3.1  (2.9) 

Plasma  adiponectin  levels, 
pg/mL 

6.0  (3.2) 

5.9  (2.8) 

6.0  (3.3) 

5.8  (3.2) 

5.7  (4.1) 

Women 

Plasma  leptin  level,  ng/mL 

<9.5 

9.6-17.7 

17.8-24.8 

24.9-37.7 

>37.8 

No.  of  controls 

148 

146 

150 

148 

149 

Age  at  blood  draw,  years 

63.3  (8.3) 

63.2  (7.9) 

64.3  (8.3) 

63.1  (7.3) 

63.1  (7.8) 

Race,  % 

White 

89.2 

95.8 

94.9 

92.1 

93.9 

Black 

2.4 

0.6 

0.9 

4.6 

4.1 

Other 

8.4 

3.6 

4.2 

3.4 

2.1 

Body-mass  index,  kg/m2 

21.8  (2.6) 

24.2  (3.9) 

25.7  (3.7) 

27.4  (3.2) 

31.3  (5.2) 

Physical  activity,  MET- 
hr/week 

21.7  (20.4) 

17.4(18.2) 

15.3  (16.8) 

13.0(13.8) 

12.0(11.8) 

Cigarette  smoking,  % 

Never 

41.1 

43.6 

52.9 

37.7 

50.7 

Past 

42.8 

44.3 

38.5 

48.7 

42.4 

20 


Current 

14.7 

10.7 

8.1 

13.0 

6.0 

Missing 

1.4 

1.5 

0.5 

0.7 

1.0 

History  of  diabetes  mellitus, 

% 

1.7 

0 

3.9 

3.9 

7.7 

Regular  multivitamin  use, 

% 

47.3 

42.6 

45.2 

43.6 

36.3 

Plasma  25(OH)D,  n mo  1/1. 

70.0  (25.3) 

65.7(21.1) 

62.7  (20.3) 

57.4  (32.4) 

54.6(18.3) 

Plasma  C-peptide  levels, 
ng/mL 

1.3  (0.6) 

1.6  (0.7) 

2.0  (1.2) 

2.2  (1.0) 

2.5  (1.2) 

Plasma  adiponectin  levels, 
pg/mL 

11.1  (6.0) 

9.6  (4.8) 

8.7  (4.7) 

7.7  (4.7) 

8.0  (4.6) 

*  Mean  (standard  deviation)  for  all  continuous  variables 

Abbreviations:  25(OH)D,  2  5 -hydroxy  vitamin  D  ;  MET-hr,  metabolic  equivalent  of  task-hour 
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Table  2.  Odds  ratios  (ORs)  and  95%  confidence  intervals  (CIs)  for  pancreatic  cancer  according  to  quintiles  of  plasma  leptin 


Quintiles  of  plasma  leptin 

_  p  * 

1  trend 

1 

2 

3 

4 

5 

Men 

Leptin  levels,  ng/mL 

Range 

<3.0 

3. 1-4.8 

4.9-7. 8 

7.9-11.7 

>11.8 

Median 

2.2 

4.0 

6.7 

9.2 

15.3 

No.  of  cases 

17 

26 

32 

27 

42 

No.  of  controls 

69 

71 

71 

72 

70 

Base  model1 

1.0 

1.57  (0.78-3.16) 

1.90  (0.95-3.80) 

1.61  (0.80-3.24) 

2.77(1.37-5.61) 

0.01 

Adjusted  model  I2 

1.0 

1.48  (0.72-3.05) 

1.80  (0.89-3.64) 

1.59  (0.78-3.25) 

2.55  (1.23-5.27) 

0.02 

Adjusted  model  II3 

1.0 

1.48  (0.72-3.05) 

1.80  (0.88-3.67) 

1.60  (0.77-3.33) 

2.54(1.13-5.72) 

0.03 

Adjusted  model  III4 

1.0 

1.68  (0.79-3.54) 

2.03  (0.96-4.29) 

1.84  (0.86-3.95) 

3.02(1.27-7.16) 

0.02 

Women 

Leptin  levels,  ng/mL 

Range 

<9.5 

9.6-17.7 

17.8-24.8 

24.9-37.7 

>37.8 

Median 

6.2 

14.0 

20.8 

30.6 

51.5 

No.  of  cases 

60 

80 

45 

60 

81 

No.  of  controls 

148 

146 

150 

148 

149 

Base  model1 

1.0 

1.36(0.91-2.03) 

0.76  (0.49-1.19) 

0.98  (0.64-1.50) 

1.27  (0.84-1.91) 

0.64 

Adjusted  model  I2 

1.0 

1.36(0.91-2.03) 

0.74  (0.48-1.16) 

0.95  (0.61-1.47) 

1.19(0.78-1.82) 

0.90 

Adjusted  model  II3 

1.0 

1.27  (0.84-1.92) 

0.68  (0.43-1.09) 

0.82  (0.50-1.33) 

0.94  (0.55-1.62) 

0.37 

Adjusted  model  III4 

1.0 

1.26  (0.82-1.93) 

0.63  (0.38-1.03) 

0.72  (0.43-1.21) 

0.84  (0.46-1.51) 

0.21 

22 


*P trend  values  were  calculated  by  the  Wald  test  of  a  score  variable  that  contained  median  values  of  quintiles. 


'ORs  and  95%  Cl  were  estimated  by  conditional  logistic  regression  conditioned  on  the  matching  factors  including  year  of  birth,  prospective  cohort 
(HPFS,  NHS,  PHS,  WHI,  WHS),  smoking  status  (never,  past,  current),  fasting  status  (fasting,  non- fasting),  and  month/year  of  blood  draw. 

2  Further  adjusted  for  race  (White,  Black,  other),  history  of  diabetes  mellitus  (yes,  no),  current  multivitamin  use  (yes,  no),  and  plasma  25(OH)D 
(continuous). 

3Further  adjusted  for  BMI  (continuous)  and  physical  activity  (continuous). 

4Further  adjusted  for  plasma  C-peptide  (continuous)  and  plasma  adiponectin  (quartiles). 

Abbreviations:  25(OH)D,  25-hydroxy  vitamin  D  ;  BMI,  body  mass  index;  HPFS,  Health  Professionals  Follow-up  Study;  NHS,  Nurses’  Health 
Study;  PHS,  Physicians’  Health  Study;  WHI,  Women’s  Health  Initiative;  WHS,  Women’s  Health  Study 
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Table  3.  Association  between  SNPs  in  the  leptin  receptor  ( LEPR )  gene  and  risk  of  pancreatic  cancer  among 
women 


SNP 

Minor 

allele 

Controls 

Cases 

Additive  Model 

N 

MAF 

(%) 

N 

MAF 

(%) 

ORa 

(95%  Cl) 

RawP- 

value 

Corrected 

P-value 

rsl0493380 

C 

646 

0.16 

272 

0.21 

1.54  (1.18-2.02) 

0.001 

0.032 

rs3790424 

G 

655 

0.31 

280 

0.24 

0.69  (0.55-0.89) 

0.003 

0.069 

rs6673324 

G 

650 

0.48 

277 

0.55 

1.38(1.11-1.71) 

0.003 

0.072 

rs2154381 

G 

653 

0.34 

281 

0.27 

0.72  (0.57-0.90) 

0.004 

0.088 

rs6662904 

A 

648 

0.42 

278 

0.48 

1.32  (1.07-1.62) 

0.009 

0.170 

rs9436746 

A 

653 

0.45 

280 

0.39 

0.75  (0.61-0.93) 

0.010 

0.180 

rs9436747 

A 

653 

0.42 

278 

0.35 

0.75  (0.60-0.93) 

0.010 

0.184 

rsll801408 

A 

657 

0.18 

283 

0.23 

1.36  (1.06-1.75) 

0.014 

0.243 

rs4 1459646 

C 

654 

0.16 

279 

0.20 

1.36  (1.05-1.77) 

0.021 

0.329 

rs2767485 

G 

650 

0.18 

280 

0.22 

1.27  (0.99-1.63) 

0.059 

0.647 

rs9436301 

G 

650 

0.23 

280 

0.27 

1.25  (0.99-1.58) 

0.063 

0.672 

rs7524834 

G 

659 

0.42 

282 

0.46 

1.21  (0.98-1.49) 

0.074 

0.724 

rsl2025906 

G 

657 

0.21 

280 

0.18 

0.79  (0.61-1.03) 

0.084 

0.770 

rs9436748 

A 

654 

0.39 

278 

0.43 

1.18  (0.95-1.46) 

0.134 

0.901 

rsl0128072 

C 

650 

0.15 

277 

0.18 

1.19  (0.90-1.58) 

0.210 

0.975 

rs7602 

A 

655 

0.20 

279 

0.23 

1.17(0.91-1.50) 

0.222 

0.981 

rsl887285 

G 

657 

0.09 

277 

0.11 

1.23  (0.88-1.72) 

0.233 

0.984 

rs3828033 

A 

647 

0.37 

276 

0.35 

0.89  (0.72-1.10) 

0.293 

0.996 

rs913199 

A 

652 

0.45 

279 

0.49 

1.09  (0.89-1.33) 

0.420 

1.000 

rsl892534 

A 

653 

0.39 

278 

0.41 

1.08  (0.87-1.34) 

0.474 

1.000 

rs379043 1 

G 

651 

0.20 

279 

0.19 

0.91  (0.71-1.18) 

0.482 

1.000 

rs3806318 

G 

654 

0.26 

279 

0.28 

1.08  (0.86-1.36) 

0.519 

1.000 

rs2148683 

G 

655 

0.47 

282 

0.49 

1.06  (0.87-1.30) 

0.537 

1.000 

rs2148682 

G 

655 

0.35 

280 

0.37 

1.05  (0.85-1.30) 

0.649 

1.000 

rsl2753193 

G 

659 

0.39 

279 

0.40 

1.05  (0.85-1.30) 

0.663 

1.000 

rs4420065 

A 

651 

0.39 

279 

0.40 

1.04  (0.84-1.29) 

0.691 

1.000 

rs4655537 

A 

654 

0.36 

275 

0.37 

1.04  (0.84-1.29) 

0.736 

1.000 

rs6700896 

A 

654 

0.39 

280 

0.40 

1.03  (0.83-1.27) 

0.776 

1.000 

rsll585329 

A 

655 

0.16 

279 

0.14 

0.97  (0.72-1.29) 

0.817 

1.000 

rs9436737 

G 

650 

0.14 

274 

0.14 

1.02  (0.76-1.36) 

0.900 

1.000 

rs3790436 

C 

653 

0.45 

277 

0.46 

1.01  (0.82-1.25) 

0.920 

1.000 

rs  1 7 12760 1 

G 

649 

0.13 

278 

0.13 

1.01  (0.75-1.37) 

0.926 

1.000 

aORs  and  95%  Cl  were  estimated  using  conditional  logistic  regression,  conditioning  on  matching  factors  including  year  of 
birth,  prospective  cohort  (HPFS,  NHS,  PHS,  WHI,  WHS)  which  also  conditions  on  gender,  smoking  status  (never,  past, 
current),  fasting  status  (fasting,  non  fasting),  and  month/year  of  blood  draw,  and  adjusted  for  race  (White,  Black,  other), 
history  of  diabetes  mellitus  (yes,  no),  current  multivitamin  use  (yes,  no),  plasma  25(OH)D  (continuous),  BMI 
(continuous),  physical  activity  (continuous),  plasma  C-peptide  (continuous)  and  plasma  adiponectin  (quartiles). 
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Abbreviations:  25(OH)D,  25-hydroxy  vitamin  D  ;  BMI,  body  mass  index;  Cl,  confidence  ratio  ;  HPFS,  Health 
Professionals  Follow-up  Study;  OR,  odds  ratio;  MAF,  minor  allele  frequency;  NHS,  Nurses’  Health  Study;  PHS, 
Physicians’  Health  Study;  SNP,  single  nucleotide  polymorphism;  WHI,  Women’s  Health  Initiative;  WHS,  Women’s 
Health  Study 
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Figure  titles  and  legends 
Figure  1. 

Title: 

Correlation  between  body-mass  index  and  plasma  leptin  in  men  and  women 

Legend: 

Scatterplot  of  body-mass  index  vs.  plasma  leptin  for  women  (blue  squares)  and  men  (red  squares).  Trendline  (line  of  best 
fit)  is  shown  as  full  (men)  or  dashed  line  (women). 

Abbreviations:  BMI,  body-mass  index  ;  P-corr,  P  correlation 

Figure  2. 

Title: 

Cohort-specific  and  meta-analysis  of  pooled  odds  ratios  for  pancreatic  cancer  according  to  plasma  leptin  levels 
(per  5  ng/mL  increase  in  plasma  leptin). 

Legend: 

Cohort-specific  multivariate  odds  ratios  conditioned  on  matching  factors  including  year  of  birth,  prospective  cohort 
(HPFS,  NHS,  PHS,  WHI,  WHS),  smoking  status  (never,  past,  current),  fasting  status  (fasting,  non-fasting),  and 
month/year  of  blood  draw,  and  adjusted  for  covariates  including  race  (White,  Black,  other),  history  of  diabetes  mellitus 
(yes,  no),  current  multivitamin  use  (yes,  no),  plasma  25(OH)D  (continuous),  body-mass  index  (continuous),  physical 
activity  (continuous),  plasma  C-peptide  (continuous),  and  plasma  adiponectin  (quartiles).  The  pooled  odds  ratio  is 
calculated  by  the  DerSimonian  and  Laird  random- effects  model.  The  solid  squares  and  horizontal  lines  correspond  to  the 
cohort-specific  multivariate  odds  ratios  and  95%  confidence  intervals,  respectively.  The  area  of  the  solid  square  reflects 
the  cohort-specific  weight  (inverse  of  the  variance).  The  open  diamond  represents  the  pooled  multivariate  odds  ratio  and 
95%  confidence  interval.  The  solid  vertical  line  indicates  an  odds  ratio  of  1.0. 

Abbreviations:  25(OH)D,  25-hydroxy  vitamin  D  ;  Cl,  confidence  ratio  ;  HPFS,  Health  Professionals  Follow-up  Study;  OR, 
odds  ratio;  NHS,  Nurses’  Health  Study;  PHS,  Physicians’  Health  Study;  WHI,  Women’s  Health  Initiative;  WHS, 
Women’s  Health  Study 
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Figure  1 
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TUMOR  METABOLISM 

Tissue  of  origin  dictates 
branched-chain  amino  acid  metabolism 
in  mutant  Kras-driven  cancers 

Jared  R.  Mayers,1’2*  Margaret  E.  Torrence,1’2*  Laura  V.  Danai,1 
Thales  Papagiannakopoulos, 1  j  Shawn  M.  Davidson,1’2  Matthew  R.  Bauer,1 
Allison  N.  Lau,1  Brian  W.  Ji,3  Purushottam  D.  Dixit,3  Aaron  M.  Hosios,1’2 
Alexander  Muir,1  Christopher  R.  Chin,1  Elizaveta  Freinkman,1’2’4’5’6  Tyler  Jacks,1’2’6 
Brian  M.  Wolpin,7  Dennis  Vitkup,3  Matthew  G.  Vander  Ileiden1’2’5’7^ 

Tumor  genetics  guides  patient  selection  for  many  new  therapies,  and  cell  culture  studies 
have  demonstrated  that  specific  mutations  can  promote  metabolic  phenotypes.  However, 
whether  tissue  context  defines  cancer  dependence  on  specific  metabolic  pathways  is 
unknown.  Kras  activation  and  Trp53  deletion  in  the  pancreas  or  the  lung  result  in  pancreatic 
ductal  adenocarinoma  (PDAC)  or  non-small  cell  lung  carcinoma  (NSCLC),  respectively,  but 
despite  the  same  initiating  events,  these  tumors  use  branched-chain  amino  acids  (BCAAs) 
differently.  NSCLC  tumors  incorporate  free  BCAAs  into  tissue  protein  and  use  BCAAs  as  a 
nitrogen  source,  whereas  PDAC  tumors  have  decreased  BCAA  uptake.  These  differences  are 
reflected  in  expression  levels  of  BCAA  catabolic  enzymes  in  both  mice  and  humans.  Loss 
of  Bcatl  and  Bcat2,  the  enzymes  responsible  for  BCAA  use,  impairs  NSCLC  tumor  formation, 
but  these  enzymes  are  not  required  for  PDAC  tumor  formation,  arguing  that  tissue  of  origin 
is  an  important  determinant  of  how  cancers  satisfy  their  metabolic  requirements. 


The  development  of  new  cancer  therapeutics 
relies  on  underlying  genetic  features  to 
identify  sensitive  patients  (J).  Mutations 
in  both  KRAS  and  TP53  are  common  genetic 
events  found  in  tumors  arising  from  many 
tissues,  and  cancers  with  these  mutations  are 
often  difficult  to  treat  (2, 3).  These  genetic  events, 
as  well  as  others  associated  with  cancer,  contrib¬ 
ute  to  the  metabolic  changes  that  support  biomass 
accumulation  and  cancer  cell  proliferation  (4). 
Oncogenic  RAS  signaling  increases  glucose  and 
glutamine  consumption  to  support  anabolic  pro¬ 
cesses  including  nucleotide,  lipid,  and  nonessential 
amino  acid  biosynthesis  and  can  also  drive  extra¬ 
cellular  protein  and  lipid  scavenging  (5).  TP53 
mutations  increase  glucose  consumption  and  gly¬ 
colytic  flux,  whereas  inactivation  of  TP53  renders 
cancer  cells  more  dependent  on  serine  uptake  and 
metabolism  (6). 

KRAS  and  TP33  mutations  are  found  in  most 
human  pancreatic  tumors  (7)  and  are  also  com¬ 
mon  in  lung  adenocarcinoma  (8).  What  is  known 
of  how  mutant  KRAS  or  disruption  of  TP53  affects 
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cancer  metabolism  is  based  on  cell  culture  studies 
in  defined  medium,  although  in  vivo  nutrient 
availability  varies  widely  between  tissues,  and 
vasculature  changes  can  limit  nutrient  access 
within  tumors  (9,  10).  The  inability  to  model 
these  differences  in  culture  has  therefore  limited 
understanding  of  how  tissue  of  origin  influences 
tumor  metabolism  (11).  Furthermore,  environ¬ 
ment  can  influence  metabolic  phenotypes  in  vitro 
(12-14),  and  metabolic  dependencies  in  vivo  can 
differ  from  those  found  in  vitro  (IS).  Metabolic 
differences  between  tumor  types  may  also  result 
from  cell-autonomous  effects,  and  tumor  meta¬ 
bolic  gene  expression  more  closely  resembles  that 
of  its  tissue  of  origin  than  that  of  other  tumors 
(16).  The  same  oncogenic  driver  can  also  cause 
different  metabolic  phenotypes  in  lung  and  liver 
tumors  (17).  This  raises  the  possibility  that  tumor 
type  is  a  major  determinant  of  some  tumor  meta¬ 
bolic  dependencies  in  vivo. 

Elevated  plasma  branched-chain  amino  acid 
(BCAA)  levels  are  found  in  early  pancreatic  ductal 
adenocarinoma  (PDAC)  but  not  in  non-small  cell 
lung  carcinoma  (NSCLC),  even  when  the  tumors 
are  initiated  by  the  same  genetic  events  (18).  To 
confirm  that  tumor  tissue  of  origin  influences 
whole-body  BCAA  metabolism,  we  used  LSL— 
KrasGI2D/*;Trp53f,°*:/JIoj!  (KP)  mice.  We  crossed 
KP  mice  to  mice  harboring  a  Cre-recombinase 
allele  driven  by  a  Pdx-1  promoter  (KP“/_C  model) 
(19)  or  delivered  viral  Cre  to  the  lungs  of  these 
mice  (20)  in  order  to  generate  models  of  PDAC 
and  NSCLC,  respectively.  Consistent  with  prior 
reports  (18),  mice  with  early  PDAC  have  increased 
levels  of  plasma  BCAAs,  whereas  mice  with  early 
NSCLC  exhibit  decreased  plasma  BCAA  levels  (fig. 
SI,  A  to  D).  When  cells  derived  from  these  tumors 
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are  implanted  subcutaneously  into  syngeneic 
hosts,  tumors  derived  from  PDAC  cells  did  not 
affect  plasma  BCAA  levels  (fig.  S1E)  (18),  whereas 
tumors  derived  from  NSCLC  cells  led  to  de¬ 
creased  plasma  BCAAs  (fig.  S1F).  These  results 
suggest  that  tumor  formation  from  NSCLC  cells 
can  cause  depletion  of  circulating  BCAAs. 

To  trace  tissue-specific  differences  in  BCAA 
metabolism  in  animals  with  pancreatic  or  lung 
tumors,  mice  were  fed  an  amino  acid-defined 
diet  in  which  20%  of  leucine  and  valine  were 
13C-labeled.  All  groups  of  mice  exhibited  similar 
levels  of  plasma  ^C-BCAA  enrichment  after  1  week 
of  exposure  to  labeled  diets  (figs.  S2,  A  and  B). 
Whereas  PDAC  tumors  contained  slightly  de¬ 
creased  free  BCAAs  relative  to  normal  pancreas, 
NSCLC  tumors  displayed  a  significant  increase 
in  labeled  free  BCAAs  compared  with  normal  lung 
(Fig.  1A  and  fig.  S2,  C  and  D).  These  differences 
are  not  a  reflection  of  different  amino  acid  com¬ 
positions  of  normal  or  tumor  tissue  in  either  the 
PDAC  or  NSCLC  models  (fig.  S3).  Because  BCAAs 
are  essential  amino  acids  that  animals  cannot 
synthesize  de  novo  (21),  these  results  suggest  that 
unlike  PDAC  tumors,  NSCLC  tumors  display  en¬ 
hanced  BCAA  uptake. 

BCAAs  have  several  potential  metabolic  fates 
in  tissues  (Fig.  IB).  They  can  be  directly  incorpo¬ 
rated  into  protein  or  reversibly  transaminated  by 
BCAA  transaminase  (Beat)  to  produce  branched- 
chain  a-ketoacids  (BCKAs)  and  glutamate.  BCKAs 
can  regenerate  BCAAs,  be  secreted,  or  be  oxida¬ 
tively  decarboxylated  by  the  branched-chain  keto- 
acid  dehydrogenase  (Bckdh)  complex  to  allow 
further  oxidation  of  the  carbon  skeleton  (21).  In 
agreement  with  increased  BCAA  uptake  in  NSCLC 
tumors,  lung  tumors  displayed  increased  labeled 
BCAA  incorporation  into  protein  compared  with 
normal  lung,  whereas  PDAC  tumors  incorporated 
less  labeled  BCAAs  relative  to  normal  pancreas 
(Fig.  1C  and  fig.  S2,  E  and  F).  Analysis  of  metabolites 
derived  from  BCAA  catabolism  revealed  that  NSCLC 
tumors  also  had  more  labeled  a-ketoisocaproate 
(KIC),  the  leucine-derived  BCKA,  whereas  no  change 
was  observed  in  levels  of  this  labeled  metabolite 
in  PDAC  tumors  (Fig.  ID  and  fig.  S2G).  No  other 
differences  in  labeled  BCAA  catabolite  levels  were 
observed  in  NSCLC  compared  with  normal  tissues, 
but  PDAC  tumors  showed  decreased  labeling  of 
the  tricarboxylic  acid  (TCA)  cycle  intermediate 
citrate  relative  to  normal  pancreas  from  labeled 
BCAAs  (Fig.  IE  and  fig.  S2,  G  to  I).  This  is  con¬ 
sistent  with  recent  work  demonstrating  minimal 
catabolism  of  BCAAs  to  TCA  intermediates  in  pro¬ 
liferating  cells  (22).  We  then  explored  whether  ex¬ 
cess  KIC  may  be  excreted  by  NSCLC  tumors  for 
further  metabolism  by  other  tissues,  such  as  liver, 
which  has  limited  Beat  but  high  Bckdh  activity 
(21, 23).  Consistent  with  this  hypothesis,  we  observed 
increased  labeling  of  downstream  leucine  metab¬ 
olites  in  the  livers  of  mice  with  lung  tumois  (fig.  S4). 
Taken  together,  these  data  suggest  that  BCAA  up¬ 
take  and  transamination,  but  not  their  subsequent 
catabolism,  may  provide  a  benefit  to  NSCLC  tu¬ 
mors,  potentially  by  acting  as  a  source  of  nitrogen. 

To  examine  whether  NSCLC  tumors,  but  not 
PDAC  tumors,  use  BCAAs  as  a  source  of  nitrogen, 
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Fig.  1.  Mice  with  NSCLC  display  increased  BCAA  uptake  and  metabolism.  (A  and  C  to  E)  Mice  were 
fed  13C-BCAA-containing  diet  for  7  days.  (A)  Relative  ion  counts  by  means  of  liquid  chromatography- 
mass  spectrometry  (LC-MS)  analysis  of  fully  labeled,  free  BCAAs  in  tumors  from  PDAC  and  NSCLC 
mice  and  normal  tissues  from  their  respective  control  mice.  Data  are  presented  as  mean  ±  SEM;  n  =  4 
control  and  n  =  4  PDAC;  n  =  4  control  and  n  =  4  NSCLC.  (B)  Diagram  of  the  leucine  catabolic  pathway. 
Red  labels  indicate  metabolites  measured  with  mass  spectrometry.  Blue  circles  indicate  13C-labeled 
carbons.  KIC,  a-ketoisocaproate.  (C)  Relative  ion  counts  by  means  of  gas  chromatography  (GC)-MS 
analysis  of  fully  labeled  BCAAs  from  protein  acid  hydrolysates  of  tumors  from  PDAC  and  NSCLC  mice 
and  normal  tissues  from  their  respective  control  mice.  Data  are  presented  as  mean  ±  SEM;  n  =  4  control 
and  n  =  4  PDAC;  n  =  4  control  and  n  =  4  NSCLC.  (D)  Relative  ion  counts  by  means  of  LC-MS  analysis  of 
fully  labeled  KIC  in  tumors  from  PDAC  and  NSCLC  mice  and  normal  tissues  from  their  respective  control 
mice.  Data  are  presented  as  mean  ±  SEM;  n  =  4  control  and  n  =  4  PDAC;  n  =  4  control  and  n  =  4  NSCLC. 
(E)  Citrate  M+2  labeling  (%)  from  [U-13C]-leucine  by  means  of  GC-MS  analysis  in  tumors  from  PDAC 
and  NSCLC  mice  and  normal  tissues  from  their  respective  control  mice.  Data  are  presented  as  mean  ± 
SEM;  n  =  4  control  and  n  =  4  PDAC;  n  =  4  control  and  n  =  4  NSCLC.  Two-tailed  t  test  was  used  for  all 
comparisons  between  two  groups.  *P  <  0.05,  **P  <  0.01,  ***P  <  0.001. 


we  fed  mice  a  modified  amino  acid  diet  in  which 
50%  of  leucine  was  labeled  with  15N,  allowing 
the  fate  of  leucine-derived  nitrogen  to  be  traced 
(Fig.  2A).  In  agreement  with  13C-tracing,  mice 
with  PDAC  demonstrated  no  differences  in  free 
^N-labeled  leucine  in  tumors  compared  with  con¬ 
trol  pancreas  (fig.  S5A)  and  had  less  ^  incorpo¬ 
ration  into  other  amino  acids  (fig.  S5B).  In  contrast, 
increased  levels  of  15N-leucine  were  found  in 


NSCLC  tumors  compared  with  normal  lung  (fig. 
S5C)  with  decreased  plasma  enrichment  of  ^N- 
leucine  in  mice  with  NSCLC  tumors  (Fig.  2B  and 
fig.  S5D).  A  relative  increase  in  15N-labeling  of 
nonessential  amino  acids,  as  well  as  valine  and 
isoleucine,  was  observed  in  both  the  free  and 
tissue-protein  amino  acid  pools  of  NSCLC  com¬ 
pared  with  control  lung  (Fig.  2C  and  fig.  S5,  C 
and  E).  Given  the  reduced  plasma  enrichment  of 
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15N-amino  acid  species  (Fig.  2B  and  fig.  S5D),  this  ^N- 
labeling  pattern  argues  that  BCAA-transamination 
mediated  by  Beat  isoforms  is  active  in  NSCLC 
tumor  tissue.  Evidence  for  increased  BCAA  trans¬ 
amination  in  NSCLC  compared  with  PDAC  cells 


-4  0  4 

is  also  evident  in  vitro  across  a  range  of  glutamine 
concentrations;  however,  tissue  culture  does  not 
recapitulate  the  same  phenotypes  observed  in 
tumors  (fig.  S6).  Downstream  of  nonessential 
amino  acid  biosynthesis,  this  nitrogen  can  also 


Fig.  2.  BCAA-derived  nitrogen  supports  non- 
essential  amino  acid  and  DNA  synthesis  in  NSCLC 
tumors.  (A)  Diagram  of  leucine  transamination  by 
Beat  and  nitrogen  (green  circles)  fate  after  transam¬ 
ination.  (B  to  D)  NSCLC  mice  were  fed  15N-leucine 
containing  diet  for  6  days.  (B)  Relative  ion  counts  by 
means  of  GC-MS  analysis  of  M+l-labeled  amino  acids 
in  plasma  of  control  and  NSCLC  mice.  Data  are  pres¬ 
ented  as  mean  ±  SEM;  n  =  5  control  and  n  =  6 
NSCLC.  (C)  Relative  ion  counts  by  means  of  GC-MS 
analysis  of  M+l-labeled  amino  acids  from  protein 
acid  hydrolysates  of  control  mouse  lung  tissue 
and  NSCLC  mouse  tumors.  Data  are  presented  as 
mean  ±  SEM;  n  =  6  control  and  n  =  6  NSCLC.  (D) 
M+1  labeling  (%)  from  15N-leucine  of  deoxynucleic 
acids  from  nucleic  acid  digest  of  control  mouse 
lung  tissue  and  NSCLC  mouse  tumors.  Data  are 
presented  as  mean  ±  SEM;  n  =  6  control  and  n  =  6 
NSCLC.  Two-tailed  t  test  was  used  for  all  compar¬ 
isons  between  two  groups.  *P  <  0.05,  **P  <  0.01, 
***P  <  0.001. 


Fig.  3.  Gene  expression  in  both  mouse  and 
human  tumors  reflects  tumor  tissue-specific 
BCAA  metabolism.  (A)  Relative  expression  of 
BCAA  metabolic  pathway  genes  in  normal  lung  and 
NSCLC  tumors  from  KP  mice.  Data  are  presented 
as  mean  ±  SEM;  n  =  6  control  and  n  =  6  NSCLC. 
(B)  Relative  expression  of  BCAA  metabolic  pathway 
genes  in  normal  pancreas  and  PDAC  tumors  from 
KP  mice.  Data  are  presented  as  mean  ±  SEM;  n  =  7 
control  and  n  =  5  PDAC.  (C)  Immunoblots  of  pro¬ 
teins  involved  in  BCAA  metabolism  in  representative 
normal  lung  and  NSCLC  tumors  (left)  and  repre¬ 
sentative  normal  pancreas  and  PDAC  tumors  (right) 
from  KP  mice.  (D)  Quantification  of  (C).  Data  are 
presented  as  mean  ±  SEM;  n  =  6  control  and  n  =  6 
NSCLC;  n  =  4  control  and  n  =  4  PDAC.  (E)  Com¬ 
parison  of  BCAA  metabolic  pathway  gene  expres¬ 
sion  in  human  NSCLC  and  PDAC  tumors  relative  to 
their  adjacent  paired  normal  tissues.  Overall  expres¬ 
sion  of  BCAA  metabolism  genes  is  significantly  de¬ 
creased  in  PDAC  (P  <  0.0001).  Two-tailed  t  test  was 
used  for  all  comparisons  between  two  groups  unless 
otherwise  stated.  *P  <0.05,  **P  <  0.01,  ***P  <  0.001. 


be  used  to  generate  nucleotides,  primarily  if  as¬ 
partate  is  synthesized  de  novo  in  these  tumors. 
Consistent  with  this  possibility,  we  found  increased 
incorporation  of  15N-label  in  both  aspartate  and 
nucleotides  (Fig.  2,  C  and  D,  and  fig.  S5E).  In 
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some  contexts,  aspartate  production  is  limiting  for 
nucleotide  biosynthesis  and  proliferation  (24, 25), 
indicating  that  BCAA  metabolism  may  be  impor¬ 
tant  for  tumor  growth. 

To  test  whether  gene  expression  differences 
might  contribute  to  differential  BCAA  metabo¬ 
lism,  we  used  quantitative  reverse  transcription 
polymerase  chain  reaction  (RT-PCR)  to  ana¬ 
lyze  mRNA  levels  in  NSCLC  and  PDAC  tumors 
compared  with  their  respective  normal  tissues. 
Consistent  with  increased  BCAA  uptake  and 
KIC  generation  in  NSCLC  tumors,  these  tumors 
displayed  increased  expression  of  the  primary 
BCAA  transporter  Slc7a5  (also  called  the  neutral 
amino  acid  transporter  Latl)  and  increased  levels 
of  Bcat2  and  Bckdh  (Fig.  3,  A,  C,  and  D).  In 
contrast,  PDAC  exhibited  decreased  expression  of 
these  genes  relative  to  normal  pancreas  (Fig.  3,  B  to 
D).  We  also  observed  increased  inhibitory  phos¬ 
phorylation  of  the  Bckdh  complex  in  lung  tumors 
(Fig.  3,  C  and  D).  Beat  expression  enables  use  of 
BCAAs  as  a  source  of  nitrogen  by  lung  tumors, 
and  inhibition  of  Bckdh  prevents  further  catab¬ 
olism  of  these  amino  acids. 

The  expression  changes  observed  in  PDAC  are 
not  specific  to  this  model;  the  related  KPC  mouse 
model  (26),  which  is  initiated  by  a  point  mutation 
in  Trp53,  showed  similar  changes  in  gene  expres¬ 
sion  (fig.  S7A).  Furthermore,  these  decreases  in 
gene  expression  do  not  appear  to  be  a  consequence 
of  the  relative  decrease  in  cancer  cellularity  of 
PDAC  tumors  (7)  because  sorted  pancreatic  can¬ 
cer  cells  showed  similar  expression  of  genes  in¬ 
volved  in  proximal  BCAA  catabolism  relative  to 
whole  tumor  extracts  (fig.  S7B).  In  further  agree¬ 
ment  with  neither  lung  nor  pancreatic  cancers 
showing  evidence  of  downstream  BCAA-carbon 
oxidation,  the  expression  of  enzymes  from  this 
pathway  was  not  markedly  different  in  either 
of  these  cancers  (fig.  S7,  C  and  D).  In  contrast, 
glycolytic  gene  expression  was  increased  in  both 
tumor  types  (fig.  S7,  E  and  F),  which  is  consistent 
with  known  increases  in  glycolysis  in  each  tumor 
type  (27-29).  Last,  to  relate  these  data  to  tissue  of 
origin,  we  performed  principal  component  and 
clustering  analyses,  which  demonstrated  segre¬ 
gation  of  each  tumor  with  the  normal  tissue  from 
which  it  arose  (figs.  S7,  G  and  H). 

To  ascertain  whether  similar  changes  in  gene 
expression  were  also  found  in  human  cancers,  we 
examined  expression  of  BCAA  catabolic  enzymes 
in  NSCLC  and  PDAC  relative  to  their  respective 
normal  tissues  in  publicly  available  data  sets  (30). 
Consistent  with  our  observations  in  mice,  human 
NSCLC  had  increased  expression  of  SLC7A5, 
BOAT,  and  BCKDH,  whereas  expression  of  BCAA 
catabolism  pathway  enzymes  was  decreased  in 
human  PDAC  (P  <  0.0001  for  the  pathway)  (Fig. 
3E  and  tables  SI  and  S2).  The  distinct  expres¬ 
sion  patterns  for  each  tumor  type  were  highly 
correlated  across  multiple  data  sets  (fig.  S7I  and 
tables  S3  to  S6).  The  similarity  between  human 
NSCLC  and  the  mouse  model  of  NSCLC  was 
observed  despite  KRAS  and  TP53  mutations 
occurring  in  <50%  of  human  tumors  (8),  and 
similar  expression  patterns  were  also  seen  in 
squamous  cell  lung  cancer  (fig.  S7I  and  table  S6), 


further  supporting  the  notion  that  tissue  of  origin 
can  dictate  metabolic  phenotype. 

The  increased  contribution  of  plasma  BCAAs 
to  biomass  in  NSCLC  tumors  suggests  that  these 
tumors  may  rely  on  BCAA  metabolism  for  growth. 
To  test  this  possibility,  we  used  clustered  regularly 
interspaced  short  palindromic  repeats  (CRISPR)- 


NSCLC  in  vitro 


Control  Clone  A  Clone  B 
Cell  Line 


NSCLC  Subcutaneous  Allografts 


NSCLC  Lung  Orthotopic  Allografts 


Control  Clone  A 
Cell  Line 


Cas9-mediated  genome  editing  to  disrupt  exon 
sequences  present  in  both  the  Bcatl  (cytosolic) 
and  the  Bcat2  (mitochondrial)  isoforms  (fig.  S8A) 
in  cancer  cell  lines  derived  from  KP  mice  with 
NSCLC  (Beat  null  clones  A  and  B)  or  PDAC  (Beat 
null)  (fig.  S8B).  Expression  analysis  and  ^-leucine 
tracing  studies  confirmed  functional  deletion  of 
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Fig.  4.  Beat  activity  is  required  for  NSCLC  tumor  growth.  (A)  Doubling  time  of  parental,  control 
CRISPR-Cas9  vector  infected  (pLenti),  and  NSCLC  Beat  null  cell  lines  in  vitro.  Data  are  presented  as  mean  ± 
SEM;  n  =  3  per  group.  Representative  experiment  from  >2  repeats.  (B)  Doubling  time  of  parental,  control 
CRISPR-Cas9  vector  infected,  and  PDAC  Beat  null  cell  lines  in  vitro.  Data  are  presented  as  mean  ±  SEM;  n  = 
3  per  group.  Representative  experiment  from  &2  repeats.  (C)  Estimated  tumor  volume  (cubic  millimeters) 
of  subcutaneous  allograft  of  control  infected  and  Beat  null  syngenic  NSCLC  cell  lines  into  C57BL/6J  mice. 
Data  are  presented  as  mean  ±  SEM;  n  =  6  per  group.  Two-way  repeated  measures  analysis  of  variance 
(AN0VA)  was  used  for  comparison  between  groups.  (D)  Estimated  tumor  volume  (cubic  millimeters)  of 
subcutaneous  allograft  of  control  infected  and  Beat-null  syngenic  PDAC  cell  lines  into  C57BL/6J  mice. 
Data  are  presented  as  mean  ±  SEM;  n  =  5  pLenti  control  and  n  =  6  Beat  null. Two-way  repeated  measures 
ANOVA  was  used  for  comparison  between  groups.  (E)  Lung  orthotopic  allograft  of  control  infected  and 
Beat-null  syngenic  NSCLC  cell  lines  into  C57BL/6J  mice;  n  =  23  vector  control  andn  =  13  Beat  null  Clone  A. 
(F)  Pancreatic  orthotopic  allograft  of  control  infected  and  Beat  null  syngenic  PDAC  cell  lines  into  C57BL/6J 
mice,  n  =  8  mice  per  group.  *P  <  0.05,  **P  <  0.01,  ***P  <  0.001. 


1164  9  SEPTEMBER  2016  •  VOL  353  ISSUE  6304 


sciencemag.org  SCIENCE 


Downloaded  from  http://science.sciencemag.org/  on  November  8,  2016 


RESEARCH  \  REPORTS 


Beat  in  both  the  NSCLC  and  PDAC  cancer  cells 
(fig.  S8,  C  to  F).  Despite  loss  of  both  Beat  isoforms, 
these  cells  proliferate  at  a  rate  that  is  similar  to 
the  parental  and  vector  control  infected  cell  lines 
in  vitro  (Fig.  4,  A  and  B).  When  Scat-null  NSCLC 
cells  were  implanted  subcutaneously  in  vivo,  how¬ 
ever,  the  ability  of  these  cells  to  form  tumors  was 
impaired,  and  one  clone  failed  to  produce  tumors 
(Fig.  4C  and  fig.  S8G).  In  contrast,  Scat-null  PDAC 
cells  implanted  subcutaneously  generated  tumors 
(Fig.  4D  and  fig.  S8H).  Additionally,  orthotopic 
transplantation  of  NSCLC  Beat- null  cells  failed 
to  form  lung  tumors  (Fig.  4E),  whereas  PDAC 
.Beat-null  cells  formed  tumors  in  the  pancreas  (Fig. 
4F).  Unlike  subcutaneously  implanted  PDAC  cells 
in  which  both  Beat  isoforms  were  knocked  out, 
these  cells  formed  smaller  tumors  in  the  pancreas 
than  did  control  cells  (fig.  S8I).  Taken  together, 
these  data  suggest  that  although  KP  lung  tu¬ 
mors  require  Beat  activity  for  growth,  this  en¬ 
zyme  activity  is  dispensable  for  KP  pancreas  tumor 
formation,  although  PDAC  tumor  growth  may  be 
aided  by  Beat  activity  in  some  tissue  environments. 

Proliferating  cells  need  to  acquire  amino  acids, 
both  to  make  protein  and  as  a  source  of  nitrogen 
for  nucleotide  and  nonessential  amino  acid  syn¬ 
thesis.  Prior  work  has  shown  that  macropinocy- 
tosis  plays  a  role  in  filling  this  requirement 
in  mutant  iMS-driven  PDAC  tumors  and  cells 
(12, 14, 31).  The  data  presented  here  argue  that 
this  process  might  be  less  active  in  mutant  Ras- 
transformed  NSCLC  tumors  that  acquire  nitrogen 
in  part  from  free  BCAAs.  Indeed,  we  observed  less 
macropinocytosis  in  cells  derived  from  mouse 
NSCLC  relative  to  mouse  PDAC  cells  (fig.  S9).  The 
decreased  reliance  of  PDAC  on  free  BCAAs,  how¬ 
ever,  does  not  necessarily  imply  that  uptake  of 
these  amino  acids  would  be  toxic  for  this  cancer. 
Overexpressing  Slc7a5  in  PDAC  cells  is  sufficient 
to  increase  leucine  uptake  (fig.  S10,  A  and  B)  but 
has  minimal  effects  on  proliferation  in  vitro  (fig. 
S10C)  or  tumor  growth  in  vivo  (fig.  S10,  D  and  E). 

A  role  for  free  BCAAs  in  supplying  nitrogen  to 
lung  cancers  is  intriguing  in  light  of  recent  studies 
in  glioblastoma  and  NSCLC  indicating  that  glu¬ 
tamine,  which  is  the  most  abundant  plasma  amino 
acid  and  serves  as  the  major  free  amino  acid 
substrate  for  nitrogen  and  carbon  in  culture  (32), 
contributes  less  to  tumor  metabolism  in  vivo 
(33, 34).  Indeed,  glucose-tracing  studies  in  humans 
and  mice  demonstrate  that  glutamine  is  net  syn¬ 
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thesized  from  glucose  (15, 33-37),  and  alternative 
sources  of  nitrogen  are  required  to  support  gluta¬ 
mine  production.  Thus,  in  these  contexts,  extrac¬ 
tion  of  nitrogen  from  BCAAs  for  de  novo  amino 
acid  and  nucleotide  biosynthesis  in  vivo  may 
explain  how  lung  tumors  satisfy  their  nitrogen 
requirements.  Consistent  with  this  possibility, 
BCAT1  expression  is  known  to  be  important  for 
glioblastoma  growth  (38),  suggesting  that  tu¬ 
mors  arising  in  tissues  other  than  the  lung  may 
also  use  BCAAs  as  a  source  of  nitrogen.  Multiple 
factors— including  local  environment,  tumor  cell 
of  origin,  and  genetic  mutations— can  lead  to 
convergent  metabolic  adaptations  in  disparate 
tumor  types. 

Elevations  in  plasma  BCAA  levels  are  associ¬ 
ated  with  early  PDAC  and  result  from  increased 
tissue  protein  breakdown  (18).  The  finding  that 
PDAC  tumors  have  decreased  use  of  circulating 
BCAAs  contributes  to  this  phenotype  as  well. 
In  contrast,  NSCLC  tumors  actively  use  BCAAs, 
leading  to  plasma  BCAA  depletions,  particularly 
because  the  liver  does  not  regulate  levels  of  these 
amino  acids  (23).  Many  patients  with  PDAC  and 
NSCLC  tumors  develop  cachexia  with  end-stage 
disease  (39).  Our  findings  suggest  that  differential 
use  of  amino  acids  by  tumors  and  the  resulting 
impact  on  whole-body  metabolism  might  play  a 
role  in  the  initiation  and  natural  history  of  cachexia. 
In  addition,  as  personalized  medicine  plays  a  larger 
role  in  the  clinical  management  of  cancer,  it  will  be 
critical  to  understand  how  cell  of  origin  and  tis¬ 
sue  environment  interact  with  genetic  events  to 
influence  metabolic  dependencies  of  tumors  and 
select  the  right  treatment  approaches  for  patients. 
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Abstract 


Background:  Many  epidemiologic  studies  are  using  metabo¬ 
lomics  to  discover  markers  of  carcinogenesis.  However,  limited 
data  are  available  on  the  influence  of  pre-analytic  blood  collection 
factors  on  metabolite  measurement. 

Methods:  We  quantified  166  metabolites  in  archived  plasma 
from  423  Health  Professionals  Follow-up  Study  and  Nurses' 
Health  Study  participants  using  liquid  chromatography-tandem 
mass  spectrometry  (LC-MS).  We  compared  multivariable-adjust¬ 
ed  geometric  mean  metabolite  LC-MS  peak  areas  across  fasting 
time,  season  of  blood  collection,  and  time  of  day  of  blood 
collection  categories. 

Results:  The  majority  of  metabolites  (160  of  166  metabo¬ 
lites)  had  geometric  mean  peak  areas  that  were  within  15% 
comparing  samples  donated  after  fasting  9  to  12  versus  >13 
hours;  greater  differences  were  observed  in  samples  donated 
after  fasting  <4  hours.  Metabolite  peak  areas  generally  were 


similar  across  season  of  blood  collection,  although  levels  of 
certain  metabolites  (e.g.,  bile  acids  and  purines/pyrimidines) 
tended  to  be  different  in  the  summer  versus  winter  months. 
After  adjusting  for  fasting  status,  geometric  mean  peak  areas  for 
bile  acids  and  vitamins,  but  not  other  metabolites,  differed  by 
time  of  day  of  blood  collection. 

Conclusion:  Fasting,  season  of  blood  collection,  and  time  of 
day  of  blood  collection  were  not  important  sources  of  variability 
in  measurements  of  most  metabolites  in  our  study.  However, 
considering  blood  collection  variables  in  the  design  or  analysis  of 
studies  may  be  important  for  certain  specific  metabolites,  partic¬ 
ularly  bile  acids,  purines/pyrimidines,  and  vitamins. 

Impact:  These  results  may  be  useful  for  investigators  formu¬ 
lating  analysis  plans  for  epidemiologic  metabolomics  studies, 
including  determining  which  metabolites  to  a  priori  exclude  from 
analyses.  Cancer  Epidemiol  Biomarkers  Prev;  25(5);  823-9.  ©2016  AACR. 


Introduction 

Metabolomic  profiling  of  archived  biospecimens  from  partici¬ 
pants  enrolled  in  epidemiologic  studies  has  the  potential  to 
facilitate  discovery  of  early  biologic  markers  of  carcinogenesis  as 
well  as  markers  of  cancer  prognosis  (1).  However,  conducting  a 
large-scale  biospecimen  collection  can  be  challenging,  particular- 
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ly  when  participants  are  spread  across  a  wide  geographic  area,  and 
standardizing  certain  aspects  of  the  blood  collection  protocol 
across  participants  may  be  particularly  difficult.  For  example,  in 
an  uncontrolled  setting,  some  participants  may  not  be  able  to 
follow  strict  protocols  with  regard  to  fasting  prior  to  blood  draw  or 
time  of  day  of  blood  collection.  Additionally,  blood  collections  in 
large  populations  may  span  several  months,  resulting  in  different 
seasons  of  blood  collection  across  participants,  and  thus  seasonal 
differences  in  exposures  (e.g.,  diet  or  light  exposure)  that  may 
affect  biomarker  levels. 

Few  studies  have  evaluated  whether  levels  of  plasma  metabo¬ 
lites  measured  by  metabolomics  profiling  platforms  vary  by  pre- 
analytic  blood  sample  collection  factors,  such  as  fasting  time,  time 
of  day  of  blood  collection,  or  season  of  blood  collection  (2-9). 
Such  information  is  necessary  to  understand  the  importance  of 
controlling  for  these  variables  in  the  design  of  metabolomics 
studies  (e.g.,  matching  cases  and  conuols  on  season  of  blood  draw 
or  restricting  the  study  population  to  those  fasting  at  least  8  hours) 
or  during  data  analysis  (e.g.,  excluding  certain  metabolites  likely 
to  be  inaccurately  measured  based  on  pre-analytic  variables). 
Failure  to  adequately  consider  these  sources  of  measurement 
error  could  lead  to  decreased  power  of  statistical  tests  and  biased 
association  measures  (10).  We  evaluated  differences  in  plasma 
metabolite  peak  areas  across  participants  in  two  large,  prospective 
cohort  studies  with  different  fasting  times  at  blood  draw,  times  of 
day  of  blood  draw,  and  seasons  of  blood  draw.  Our  analyses 
considered  data  on  166  metabolites,  measured  by  a  liquid  chro¬ 
matography-tandem  mass  spectrometry  (LC-MS)  metabolomics 
platform,  which  performed  well  in  our  previous  pilot  testing  (11). 
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Materials  and  Methods 

Study  population 

The  Nurses'  Health  Study  (NHS)  was  established  in  1976 
when  121,700  female  nurses  ages  30  to  55  years  completed  a 
mailed  questionnaire.  The  Health  Professionals  Follow-up 
Study  (HPFS)  was  initiated  in  1986  when  51,529  men  ages 
40  to  75  years  and  working  in  health  professions  completed  a 
mailed  questionnaire.  Both  cohorts  enrolled  only  adults  living 
in  the  United  States  and  have  followed  participants  since 
enrollment  via  mailed  biennial  questionnaires.  Blood  samples 
were  collected  in  heparin  tubes  from  32,826  women  in  NHS 
from  1989  to  1990  and  in  EDTA  tubes  from  18,225  men  in 
HPFS  from  1993  to  1995.  Participants  arranged  to  have  their 
blood  collected  and  shipped  the  samples  via  overnight  mail 
with  an  icepack  to  a  central  laboratory.  On  arrival,  the  samples 
were  processed  and  aliquots  of  plasma,  white  blood  cells,  and 
red  blood  cells  were  stored  in  liquid  nitrogen  freezers;  >95%  of 
samples  arrived  within  24  hours  of  collection  (12,  13).  Parti¬ 
cipants  provided  information  on  the  date  of  blood  collection, 
time  at  which  blood  was  drawn,  and  number  of  hours  since  last 
eating  before  blood  draw  on  a  questionnaire  returned  with  the 
samples. 

Our  analyses  included  423  participants  (252  from  NHS  and 
171  from  HPFS)  selected  as  controls  for  a  previous  study  of 
plasma  metabolomics  and  pancreatic  cancer  (14).  All  partici¬ 
pants  provided  implied  consent  by  returning  the  question¬ 
naires  and  blood  samples.  The  study  was  approved  by  the 
Human  Research  Committee  at  Brigham  and  Women's  Hospi¬ 
tal  (Boston,  MA). 

Metabolite  profiling 

Profiles  of  231  endogenous  polar  metabolites  and  lipids  were 
obtained  using  LC-MS  at  the  Broad  Institute  of  the  Massachusetts 
Institute  of  Technology  and  Harvard  University  (Cambridge, 
MA).  A  detailed  description  of  the  metabolite  profiling  methods 
has  been  previously  published  (14).  For  the  current  analysis,  we 
excluded  3  5  metabolites  that  were  not  reasonably  reproducible  in 
samples  processed  up  to  24  hours  after  blood  draw  (due  to  the 
design  of  our  collections)  in  previous  pilot  tests  (11,  14).  Rea¬ 
sonable  reproducibility  was  defined  as  an  intraclass  correlation 
(ICC)  >0.75  comparing  samples  processed  immediately  versus 
24  hours  later  ( 1 5 ) .  In  addition,  we  excluded  1 5  metabolites  with 
inter-assay  coefficients  of  variation  >25%,  as  measured  in  blinded 
quality  control  samples,  and  15  metabolites  with  undetectable 
levels  in  >10%  of  participants.  Thus,  166  metabolites  were 
retained  for  these  analyses. 

Statistical  analysis 

Metabolite  levels  were  reported  as  measured  LC-MS  peak  areas, 
which  are  proportional  to  metabolite  concentration.  We  grouped 
metabolites  into  9  categories:  amines  (n  =  5);  amino  acid  deri¬ 
vatives  (n  =  13);  amino  acids  (n  =  28);  bile  acids  ( n  =  5);  lipids 
and  lipid  metabolites  (n  =  88);  organic  acids  (n  =  11);  purines, 
pyrimidines,  and  derivatives  (n  =  3);  vitamins  (n  =  3);  and  other 
(n  =  10;  see  Supplementary  Table  SI  for  the  list  of  metabolites 
within  each  category) .  To  quantify  the  strength  of  the  relationship 
between  metabolites  in  each  category,  after  adjusting  for  potential 
differences  between  participants  that  might  affect  metabolite 
levels,  we  calculated  Spearman  partial  correlation  coefficients. 
The  covariates  included  age  at  blood  draw  (years,  continuous), 
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race  (white,  non-white),  body  mass  index  at  blood  draw  (con¬ 
tinuous,  kg/m2),  physical  activity  (metabolic  equivalent  hours/ 
week,  continuous),  sex/menopausal  status/hormone  use  (male, 
female/premenopausal,  female/postmenopausal  and  no  hor¬ 
mone  therapy  use,  female/postmenopausal  and  hormone  thera¬ 
py  use),  and  total  caloric  intake  (kcal,  continuous).  Data  on 
covariates  were  obtained  from  the  cohort  questionnaire  closest 
to  the  date  of  blood  collection. 

For  the  analyses  of  metabolite  levels  according  to  blood 
sample  collection  factors,  each  metabolite  peak  area  was  log- 
transformed  to  improve  the  normality  of  its  distribution  and 
included  as  the  dependent  variable  in  a  robust  variance 
linear  regression  model  (PROC  MIXED,  version  9.3;  SAS 
Institute).  The  primary  independent  variable  in  the  model 
was  the  pre-analytic  blood  collection  variable  of  interest: 
fasting  time  [<4  (referent),  5-8,  9-12,  >13  hours],  season 
of  blood  draw  [February-April  (referent),  May-July,  August- 
October,  November-January],  or  time  of  day  of  blood  draw 
[6-8  am  (referent),  9-10  am,  11  am-12  pm,  1-2  pm,  3  pm- 
12  am].  We  calculated  the  mean  log-transformed  metabolite 
peak  area  in  each  category  of  the  blood  collection  variable 
after  adjusting  for  the  covariates  (i.e.,  marginal  means).  We 
then  calculated  multivariable-adjusted  geometric  mean  peak 
areas  by  taking  the  anti-log  of  each  marginal  mean.  We 
calculated  the  percentage  difference  in  geometric  mean  metab¬ 
olite  peak  area  between  each  category  of  the  blood  collection 
variable  and  the  referent  category  (e.g.,  fasting  9-12,  5-8,  or 
<4  hours  vs.  >13  hours)  by  subtracting  the  referent  geometric 
mean  peak  area  from  the  geometric  mean  peak  area  of  the 
category  of  interest  and  dividing  by  the  referent  geometric 
mean  peak  area. 

Within  each  metabolite  category,  we  calculated  the  median  as 
well  as  the  10th  and  90th  percentiles  of  the  percentage  differences 
in  the  geomeuic  mean  metabolite  peak  area  for  each  blood 
collection  variable  category  (except  the  referent  group).  We  also 
calculated  the  percentage  (95%  confidence  interval)  of  metabo¬ 
lites  in  each  metabolite  category  whose  geometric  mean  peak 
areas  were  ±15%  of  the  referent  category  geomeuic  mean  and 
considered  alternate  cutoff  points  (10%,  20%)  in  sensitivity 
analyses. 

Because  fasting  time  and  time  of  day  of  blood  draw  are  highly 
correlated,  analyses  of  time  of  day  of  blood  draw  were  further 
adjusted  for  fasting  time  (<4,  5-8,  9-12,  >13  hours).  Also,  to 
further  explore  the  impact  of  fasting  time  on  associations  between 
time  of  day  of  blood  draw  and  metabolite  concenuation,  we 
conducted  separate  analyses  among  individuals  fasting  <8  hours 
and  those  fasting  >8  hours.  Due  to  the  smaller  sample  size,  time  of 
day  of  blood  draw  was  evaluated  in  three  categories  in  these 
analyses  (6-8  am,  9  am-12  pm,  1  pm-12  am). 

Results 

The  mean  age  at  blood  draw  was  65  years  among  men  in  HPFS 
and  60  years  among  women  in  NHS  (Table  1 ) .  Among  both  men 
and  women,  a  large  majority  were  white  (94%  in  HPFS;  99%  in 
NHS).  The  most  common  season  of  blood  collection  was  late 
summer/early  autumn  among  men  (46%  in  August-October) 
and  late  winter/early  spring  among  women  (36%  in  February- 
April).  The  majority  of  participants  in  both  cohorts  had  fasted  for 
more  than  8  hours  at  blood  collection  and  had  their  blood  drawn 
in  the  morning. 
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Pre-analytic  Blood  Collection  Factors  and  Metabolomics 


Table  1.  Age-standardized  characteristics  of  the  study  population3 


Characteristic 

HPFS 
(n  =  171) 

NHS 

(n  =  252) 

Age  at  blood  draw,  y,  mean  (SD) 

65.2  (7.8) 

60.2  (6.2) 

White,  % 

94.4 

99.4 

Height,  inches,  mean  (SD) 

70.1  (2.9) 

64.4  (2.4) 

Body  mass  index,  kg/m2,  mean  (SD) 

26.0  (3.2) 

25.9  (5.0) 

Physical  activity,  MET-hr/wk,  mean  (SD) 

38.0  (35.1) 

17.7  (20.3) 

Season  of  blood  draw,  % 

Feb-Apr 

10.0 

36.2 

May-Jul 

27.7 

31.4 

Aug-Oct 

45.9 

15.1 

Nov-Jan 

16.4 

17.3 

Fasting  time  at  blood  draw,  % 

<4  h 

31.9 

18.5 

5-8  h 

11.6 

5.0 

9-12  h 

36.8 

45.7 

>13  h 

197 

30.8 

Time  of  day  of  blood  draw,  % 

6-8  am 

29.9 

37.7 

9-10  am 

37.0 

45.7 

11  am-noon 

9.9 

7.6 

1-2  pm 

11.6 

5.2 

3  pm-midnight 

11.6 

3.8 

Abbreviations:  MET,  metabolic  equivalent;  SD,  standard  deviation. 

aAII  values  (except  age  at  blood  draw)  are  standardized  to  the  age  distribution  of 

the  study  population. 


The  distribution  of  interassay  coefficients  of  variation  for 
metabolites  in  each  category  is  shown  in  Table  2.  Spearman 
partial  correlation  coefficients  among  amines,  amino  acid  deri¬ 
vatives,  and  purines  and  pyrimidines  were  weak  to  moderate 
(Supplementary  Table  S2).  Correlations  among  amino  acids,  bile 
acids,  lipids,  and  organic  acids  ranged  from  nearly  0  to  >0.8. 

Metabolite  peak  areas  were  similar  across  categories  of  fasting 
time  for  the  majority  of  metabolites,  particularly  among  people 
who  had  last  eaten  >4  hours  prior  to  blood  draw  (Table  3). 
Specifically,  geometric  mean  metabolite  peak  areas  measured  in 
samples  collected  after  9  to  12  or  5  to  8  hours  fasting  were  within 
15%  of  geometric  mean  metabolite  peak  areas  in  samples  col¬ 
lected  after  fasting  at  least  13  hours  for  91%  to  100%  of  amino 
acid  derivatives,  amino  acids,  lipids  and  lipid  metabolites,  organic 
acids,  and  purines  and  pyrimidines.  Metabolite  peak  areas  were 
less  similarwhen  comparing  samples  from  participants  fasting  <4 
versus  >13  hours.  In  particular,  geometric  mean  peak  areas 
comparing  participants  fasting  <4  versus  >13  hours  were  within 
15%  for  only  1  of  5  bile  acids  and  none  of  the  3  vitamins. 
Supplementary  Table  SI  shows  the  percentage  differences  in 
geometric  mean  peak  areas  by  fasting  time  for  each  metabolite. 
Results  for  the  10%  and  20%  cutoffs  were  similar  to  the  15% 
cutoff  (Supplementary  Table  S3). 

Geometric  mean  peak  areas  of  most  metabolites  were  similar 
across  seasons  of  blood  collection  (Table  4  and  Supplementary 
Table  S4).  Among  participants  whose  blood  was  collected  in 
May-July,  August-October,  or  November-January,  geometric 
mean  peak  areas  for  at  least  80%  of  amines,  amino  acid 
derivatives,  amino  acids,  and  lipids  and  lipid  metabolites  were 
within  15%  of  those  among  participants  whose  blood  was 
collected  between  February  and  April.  However,  levels  of  cer¬ 
tain  types  of  metabolites  appeared  to  vary  by  season.  For 
example,  geometric  mean  metabolite  peak  areas  in  samples 
collected  during  peak  sun  months  (August-October)  were 
within  15%  of  those  collected  during  low  sun  months 


(February-April)  for  just  67%  of  purines/pyrimidines,  55%  of 
organic  acids,  and  none  of  the  bile  acids.  When  using  the  more 
conservative  criterion  of  percentage  difference  <10%,  these 
figures  were  33%,  55%,  and  0%,  respectively  (Supplementary 
Table  S5). 

We  investigated  percentage  differences  in  geometric  mean 
metabolite  peak  areas  by  time  of  day  of  blood  collection 
adjusting  for  fasting  time  (Table  5  and  Supplementary  Tables 
S6  and  S7).  Geometric  mean  levels  of  amino  acid  derivatives, 
amino  acids,  lipids  and  lipid  metabolites,  organic  acids,  and 
purines  and  pyrimidines  generally  were  similar  across  indivi¬ 
duals  whose  blood  was  collected  at  different  times  during  the 
day.  However,  levels  of  certain  metabolites,  such  as  bile  acids 
and  vitamins,  appeared  to  be  more  variable  during  the  day.  To 
further  examine  these  patterns,  we  conducted  additional  anal¬ 
yses  separately  among  participants  fasting  <8  hours  and  parti¬ 
cipants  fasting  >8  hours  (Supplementary  Table  S8).  Differences 
in  geometric  mean  metabolite  peak  areas  for  bile  acids  by  time 
of  day  of  blood  draw  were  evident  among  both  participants 
fasting  >8  and  <8  hours.  Among  participants  fasting  >8  hours, 
bile  acid  levels  tended  to  be  lower  in  those  whose  time  of  blood 
draw  was  9  am- 12  pm  or  1  pm- 12  am  versus  6-8  am  (per¬ 
centage  differences,  10th  to  90th  percentiles  =  -25%,  -32%  to 
-14%,  and  -15%,  -52%  to  -4%,  respectively).  Differences  in 
vitamin  peak  areas  by  time  of  day  of  blood  draw  were  mainly 
apparent  among  participants  fasting  <8  hours;  vitamin  levels 
tended  to  be  higher  among  those  who  blood  was  drawn 
between  9  am  and  12  am  versus  6-8  am. 

Discussion 

Overall,  our  results  suggest  that,  after  excluding  metabolites 
that  are  not  stable  with  delays  in  processing  up  to  24  hours,  fasting 
time,  time  of  day  of  blood  draw,  and  season  of  blood  draw  are  not 
important  sources  of  variability  for  the  majority  of  metabolites 
measured  by  our  metabolomics  profiling  platform,  particularly 
amines,  amino  acids,  amino  acid  derivatives,  and  lipids  and  lipid 
metabolites.  However,  our  results  also  suggest  that  careful  atten¬ 
tion  to  pre-analytic  variables  related  to  blood  collection  may  be 
necessary  in  analyses  focused  on  certain  types  of  metabolites. 
Specifically,  we  observed  variability  in  levels  of  organic  acids, 
purines  and  pyrimidines,  bile  acids,  and  vitamins  by  fasting 


Table  2.  Distribution  of  inter-assay  coefficients  of  variation  for  metabolites  in 
each  category3 


Metabolite  category 

Metabolites,  n 

Median 

(minimum  -  maximum) 
interassay  coefficient  of 
variation  (%) 

Amines 

5 

8.9  (8.2-22.5) 

Amino  acid  derivatives 

13 

13.0  (9.6-23.9) 

Amino  acids 

28 

10.4  (6.9-18.0) 

Bile  acids 

5 

17.6  (12.1-20.5) 

Lipids  and  lipid  metabolites 

88 

9.6  (4.0-25.0) 

Organic  acids 

11 

13.3  (5.2-24.8) 

Purines,  pyrimidines, 
and  derivatives 

3 

9.0  (6.8-20.5) 

Vitamins 

3 

14.8  (11.5-16.1) 

Other6 

10 

10.5  (6.6-16.3) 

aMetabolites  with  an  interassay  coefficient  of  variation  >25%  were  excluded 
from  analysis. 

bOther  metabolites  include  allantoin,  cotinine,  creatine,  creatinine,  hippurate, 
hydroxyphenylacetate,  phenylethylamine,  salicylurate,  sorbitol,  and  uracil. 
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status,  especially  when  comparing  participants  who  ate  relatively 
recently  (fasting  time  <4  hours)  to  those  whose  last  meal  was 
more  distant  (fasting  time  >13  hours).  Regarding  differences  by 
season  of  blood  draw,  we  observed  a  pattern  in  which  measure¬ 
ments  of  bile  acids,  organic  acids,  and  purines  and  pyrimidines 
were  most  variable  between  months  with  peak  sun  (May-Octo- 
ber)  versus  low  sun  (February-April)  in  the  Northern  Hemi¬ 
sphere,  where  the  large  majority  of  participants  in  our  study 
reside.  Finally,  after  adjusting  for  fasting  time,  having  blood 
drawn  in  the  late  morning,  afternoon,  or  evening  compared  with 
the  early  morning  (6-8  am)  was  associated  with  variability  in 
measured  peak  areas  for  bile  acids  and  vitamins,  although  vari¬ 
ation  in  vitamin  levels  by  time  of  day  of  blood  draw  was  restricted 
mainly  to  individuals  fasting  <8  hours. 

Few  previous  studies  have  investigated  differences  in  metab¬ 
olite  profiles  by  pre-analytic  variables  related  to  blood  collec¬ 
tion.  Brauer  and  colleagues  conducted  amino  acid  metabolite 
profiling  using  electrospray  tandem-mass  spectrometry  in 
whole  blood  collected  from  10  adults  after  fasting  >5  hours 
as  well  as  3  and  5  hours  after  eating  a  standardized  meal  (5). 
They  observed  changes  in  amino  acid  levels  3  and  5  hours 
postprandial,  suggesting  that  a  fasting  period  longer  than  5 
hours  should  be  considered  for  measurement  of  amino  acids. 
Similarly,  in  our  analysis,  fasting  time  between  5  and  12  hours 
compared  with  >13  hours  appeared  to  have  minimal  influence 
on  levels  of  amino  acids.  We  observed  more  differences  in 
amino  acid  levels  when  comparing  samples  from  participants 
fasting  <4  hours  versus  >13  hours,  although  percentage  differ¬ 
ences  in  geometric  mean  peak  areas  remained  within  15%  for 
the  large  majority  (86%)  of  amino  acids.  In  our  study,  bile  acids 
and  vitamins  were  the  most  variable  types  of  metabolites  across 
fasting  time  categories  and  tended  to  be  higher  in  participants 
fasting  <4  hours  versus  >13  hours.  This  pattern  likely  reflects 
the  intake  of  nutrients  and  reabsorption  of  bile  acids  from  the 
intestine  following  ingestion  of  food  (16).  In  general,  our 
results  regarding  variability  in  metabolite  levels  by  fasting 
status  are  congruent  with  the  conclusion  of  Sampson  and 
colleagues  that  fasting  accounts  for  little  of  the  variability  in 
metabolites  and,  therefore,  metabolomics  can  be  a  useful  tool 
even  in  studies  without  strict  guidelines  on  fasting  before  blood 
draws  (4). 

Few  data  have  been  published  on  potential  seasonal  differences 
in  metabolite  profiles.  The  differences  we  observed  in  levels  of  bile 
acids,  organic  acids,  and  purines  and  pyrimidines  across  indivi¬ 
duals  with  different  seasons  of  blood  draw  might  reflect  seasonal 
variation  in  diet  or  other  exposures,  such  as  sunlight.  Interestingly, 
the  pattern  we  observed  of  higher  urate  levels  among  participants 
who  donated  blood  in  the  summer  versus  winter  months  has  been 
previously  reported,  although  the  mechanism  underlying  this 
pattern  is  not  clear  (17).  Overall,  season  did  not  appear  to  be  an 
important  source  of  variability  for  most  types  of  metabolites  in 
our  study. 

Several  studies  have  reported  on  diurnal  variation  in  biomarker 
concentrations.  Ang  and  colleagues  observed  significant  24-hour 
variation  in  six  amino  acids  (leucine,  lysine,  methionine,  phe¬ 
nylalanine,  proline,  and  tyrosine)  in  a  laboratory  study  of  8  men 
who  consumed  a  nutritional  drink  every  hour  while  awake  (3). 
Similarly,  in  our  analysis  of  individuals  who  fasted  <8  hours 
before  blood  draw,  the  percentage  difference  in  geometric  mean 
metabolite  peak  area  comparing  individuals  who  donated  blood 
between  1pm  and  1 2  am  with  those  who  donated  blood  between 
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6  am  and  8  am  was  >15%  for  four  of  the  six  amino  adds 
highlighted  by  Ang  and  colleagues  (data  not  shown).  However, 
we  did  not  observe  large  differences  in  the  six  amino  acids  by 
time  of  day  of  blood  draw  among  participants  who  had  fasted 
more  than  8  hours  at  blood  draw  (all  percentage  differences 
<10%,  data  not  shown).  The  differences  we  observed  in  bile 
acids  by  time  of  day  of  blood  draw  among  participants  fasting 
<8  hours  are  consistent  with  known  changes  in  concentrations 
of  circulating  bile  acids  following  meal  ingestion.  In  addition, 
serum  levels  of  bile  acids  have  been  shown  to  decrease  during 
the  day  among  fasting  individuals  (18),  consistent  with  our 
results  among  individuals  fasting  >8  hours.  The  differences  we 
observed  in  levels  of  vitamins  by  time  of  day  of  blood  draw, 
only  among  participants  fasting  <8  hours,  also  likely  reflect 
food  ingestion.  However,  it  should  be  noted  that  the  three 
vitamins  included  in  our  analysis  were  B-complex  vitamins, 
which  are  water-soluble  and  not  stored  in  the  body.  Thus,  our 
results  likely  do  not  reflect  variability  in  levels  of  fat-soluble 
vitamins  during  the  day. 

Several  limitations  of  our  study  should  be  considered.  First,  we 
examined  differences  in  metabolite  levels  across  participants  with 
different  blood  collection  characteristics,  rather  than  differences 
within  individuals  observed  under  different  conditions.  However, 
we  adjusted  the  mean  metabolite  peak  areas  for  a  variety  of 
individual-level  characteristics  (e.g.,  sex,  age,  and  body  mass 
index)  that  might  influence  metabolite  levels  and  our  findings 
are  generally  consistent  with  those  observed  in  the  few  clinic- 
based  studies.  Still,  it  is  possible  that  some  of  our  novel  findings, 
particularly  those  related  to  season  of  blood  draw,  are  due  to 
chance  and  should  be  confirmed  in  future  studies.  In  addition,  we 
excluded  35  metabolites  from  analysis  that  were  not  stable  with 
delays  in  processing  up  to  24  hours,  as  demonstrated  in  our 
previous  pilot  testing  ( 1 1 ) ,  in  order  to  focus  on  metabolites  likely 
to  be  informative  in  our  study  population.  Thus,  some  of  the 
metabolite  categories  in  our  analysis  contained  small  numbers  of 
metabolites  (e.g.,  purines  and  pyrimidines,  vitamins)  and  certain 
categories  were  excluded  entirely  (e.g.,  carbohydrates);  this  limits 
our  ability  to  make  conclusions  about  the  influence  of  blood 
collection  characteristics  on  certain  metabolite  types.  Finally, 
several  factors  may  limit  the  generalizability  of  our  results.  Our 
primary  analyses  summarized  results  within  metabolite  catego¬ 
ries.  Correlations  between  metabolites  within  categories  will 
influence  the  amount  of  variability  in  percentage  differences 
observed  (Supplementary  Table  S2).  Therefore,  our  results  may 
not  be  generalizable  to  other  metabolomics  platforms  in  which 
correlations  between  metabolites  differ  from  those  in  our  study  or 
platforms  that  measured  different  sets  of  metabolites.  Also,  our 
results  may  not  apply  to  study  populations  that  differ  significantly 
from  the  generally  healthy,  middle-  and  older-aged  adults  in  our 
study. 

In  conclusion,  the  results  of  our  analysis  suggest  that  fasting 
time,  season  of  blood  draw,  and  time  of  day  of  blood  draw  are 
not  important  sources  of  variability  for  several  major  classes  of 
metabolites  measured  by  an  LC-MS  metabolomics  platform, 
including  amines,  amino  acid  derivatives,  amino  acids,  and 
lipids  and  lipid  metabolites.  Thus,  epidemiologic  studies  with 
archived  blood  specimens  may  be  useful  settings  for  future 
metabolomics  research,  even  if  strict  blood  collection  protocols 
were  not  instituted  or  data  on  blood  collection  characteristics 
are  not  available.  However,  our  results  also  suggest  that  sample 
collection  factors  may  be  important  to  consider  if  certain 


metabolite  classes  are  of  particular  scientific  interest.  Investi¬ 
gation  of  water-soluble  vitamins  should  be  restricted  to  indi¬ 
viduals  fasting  at  least  8  hours.  In  addition,  evaluating  a 
potential  influence  of  season  of  blood  draw  on  vitamin  levels 
may  be  prudent  (e.g.,  in  sensitivity  analyses).  Due  to  variability 
by  season  and  time  of  day  of  blood  draw,  analyses  of  bile  acids 
may  be  difficult  in  epidemiologic  studies.  If  bile  acids  are  of 
particular  interest,  analyses  should  be  limited  to  participants 
fasting  at  least  8  hours  at  blood  draw.  Overall,  our  results  may 
be  useful  for  epidemiologic  studies  to  a  priori  exclude  meta¬ 
bolites  subject  to  substantial  measurement  error  from  analyses, 
and  thus  reduce  penalties  for  multiple  testing,  and  when 
considering  variables  related  to  blood  collection  to  include  in 
analysis  plans. 
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